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ABSTRACT: 

Requirements for indirect computation of geostrophic surface currents 
over large ocean areas are discussed. These requirements point to a need 
to simplify standard geostrophic computations, and to separate the first 
order thermal and haline contributions to geostrophic flow. 

The equations of motion for geostrophic flow are reviewed and the 
standard geostrophic computations discussed. Errors and limitations in the 
geostrophic method are reviewed. Previous attempts to simplify geostrophic 
computations are discussed and shown to be inadequate for synoptic compu- 
tation over hemispheric ocean areas. 

It is shown that the Helland-Hansen equation can be rewritten such 
that the geostrophic velocity is composed of a temperature structure term 
and a salinity structure term. In order to apply this modified equation to 
hemispheric synoptic geostrophic computations, simple expressions are 
required for the dependence of specific volume on temperature at constant 
salinity and pressure, (Ba/BT)g p, and on salinity at constant temperature 
and pressure, (Ba/dS^p . Two approaches are explored to derive these 
expressions, using: 

1) Experimental P-V-T-S data collected in the laboratory. 

2) P-V-T-S data derived from an equation of state for sea water. 

Numerical fitting of the experimental P-V-T-S data shows that the de- 
pendence of specific volume on temperature, a(T)g p , can be expressed 
by a quadratic equation and the dependence of specific volume on salinity, 
a(S)'p / P/ can be expressed by a linear equation. The coefficients in the 
temperature dependence equation are a function of salinity and pressure, 
and the coefficient in the salinity dependence is a function of temperature 
and pressure. However, it is shown that the lack of experimental data and 
the presence of small errors in the data lead to inconsistencies in the values 
of the coefficients derived. This difficulty is corrected using Ekman's equa- 
tion of state for sea water to generate P-V-T-S data . 

Using derived expressions for ( Ba/3T)g p and (Ba/3S)„ „ an equation 
is derived from the Helland-Hansen relationship and is calleathe Tempera- 
ture-Salinity Gradient scheme for computing geostrophic currents, or simply 
T-S Gradient scheme. This scheme gives equivalent results to the standard 
geostrophic computations, yet is computationally much simpler and faster 
than the standard scheme. Other advantages are (1) the variables observed 
in the ocean, temperature and salinity structure are used directly in the com- 
putations (2) no other quantities such as Cf andcr^ need be computed, and 



(3) thermal and haline contributions to the geostrophic flow can be expressed 
independently. 



The T-S Gradient scheme is applied to existing hemispheric fields of 
T and S data. It is shown that in the expressions for (da/BT)g p and 
(Sa/dS) T p the use of fixed coefficients over the entire Pacific Ocean 
relative to a reference level of 1000 db, introduces less than five percent 
error in the surface velocity. This application shows the feasibility of 
making geostrophic computations by the T-S gradient scheme. Unfortunately, 
it also shows the deficiencies of the available fields for hemispheric compu- 
tations by any scheme. These deficiencies are (1) the inaccuracy of the 
data, (2) the large grid spacing, and (3) the inconsistency of the tempera- 
ture and salinity fields. These must be corrected before hemispheric geo- 
strophic computations can be made. The construction of new fields to meet 
the requirements of geostrophic computations using available data is dis- 
cussed . 

Two other aspects of current computation were investigated. First, 
using the capability of the T-S Gradient scheme to separate flow into thermal 
and haline components, the relationship between the thermal component and 
the total geostrophic surface current in the Gulf Stream water mass in the 
Grand Banks was studied. Results show that these quantities are satis- 
factorily related by a linear expression which allows the determination of 
the total geostrophic surface flow in this water mass using only temperature 
measurements. Such techniques may be useful in reducing survey time of 
the Coast Guard Ice Patrol in this region. 

A second aspect investigated is the relationship between the indirectly 
computed currents and measured currents in the California current. Two 
experiments were performed in which drogue-measured surface currents were 
compared to indirectly computed currents given by the Fleet Numerical Weather 
Central and those computed by the Helland-Hansen and T-S Gradient schemes 
using standard hydrographic data. This investigation points out the signifi- 
cant need for more work on this aspect. 
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THE T-S GRADIENT METHOD, A NEW METHOD 
OF COMPUTING GEOSTROP HIC CURRENTS 
OVER LARGE OCEAN AREAS 



I. INTRODUCTION 

The purpose of the research reported in this study is to improve 
the contemporary techniques of computing synoptic geostrophic ocean 
surface currents over large ocean regions. Techniques that are used 
today are either gross approximations with little scientific basis or 
applicable to only small regions, or both. Yet knowledge of surface 
currents is important, and in some cases essential to public, com- 
mercial, and military use of the sea. 

The U. S. Coast Guard uses surface current information in 
search and rescue, as well as in the forecasting of iceberg drift for 
the shipping lanes in the Northwest Atlantic (Lenczyk, 1964). The 
U. S. Navy Fleet Numerical Weather Central, Monterey, California, 
computes surface currents every 12 hours over the Northern Hemi- 
sphere principally to aid in thermal structure forecasting (Hubert 
and Laevastu, 1967). The U. S. Naval Oceanographic Office computes 
surface currents in the ASWEPS (Antisubmarine Warfare Environ- 
mental Prediction System) area of the Northwest Atlantic for thermal 
structure forecasting (James, 1966). 

Mariners have been using knowledge of ocean currents probably 



2 



since the beginning of navigation- As early as 600 BC the east to 
west ei rcum -navigation of Africa may have been aided by favorable 
winds and currents (von Arx, 1962). Early Greek and Arab traders 
could penetrate the Arabian Sea and Bay of Bengal because of 
favorable currents under the influence of the monsoon winds. Spanish 
trade between Mexico and the Philippines was also aided by knowledge 
of the winds and currents (Jones, 1939). The Spanish captains would 
take advantage of the North Equatorial Current on their westerly 
journey, and return to Mexico along a more northerly route reaching 
land fall around Cape Mendocino. An early chart of the Gulf Stream 
was published by Benjamin Franklin in 1770. This chart is said to 
have improved mail service from England to the American Colonies 
by two weeks per crossing (Groen, 1967). Pilot charts, published by 
the Naval Oceanographic Office, provide monthly summaries of wind 
and surface current conditions. Surface current atlases giving the 
monthly mean values are also available from the Naval Oceanographic 
Office. 

Flow in the surface layers of the deep ocean is influenced by 
several forces: the most important are the horizontal pressure 

gradient and viscous forces within the fluid, the action of wind stress, 
tidal forces, and the Coriolis force. Today incomplete knowledge of 
ocean dynamics and lack of suitable available data restricts the in- 
direct determination of surface currents to the computation of the 
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geostrophic and wind driven components (Hubert, 1965). 

Computation of geostrophic flow over a large ocean area on a 
synoptic basis are facilitated if the computations be as simple and as 
fast as possible. Standard geostrophic computations are unsuitable 
for many reasons: 

1) The computations are numerous and consume more computer 
time than is available for this task. 

2) Several quantities must be computed that are not used beyond 
the analysis of currents. 

3) The temperature and salinity data that are required for 
standard geostrophic computations are not yet available in the neces- 
sary accuracy and density. 

4) It would be useful to compute the temperature and salinity 
structure contributions to the geostrophic current independently. 

This is difficult by the standard method. 

Further explanation ofthe last two points is in order. Because 
of the relative ease of measuring temperature structure over the 
salinity structure, temperature is a much more commonly observed 
variable in the ocean. Until the recent development of in situ con- 
ductivity measurement, rapid determination of salinity structure was 
not possible. For this reason, while nearly synoptic thermal struc- 
ture information may be available from many areas of the oceans, 
salinity structure is available only from the historical records 
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contained in oceanographic data centers and various atlases derived 
trom this historical data. It will probably be several years before 
salinity structure is as commonly observed as thermal structure. 

It is because of this compatability that it is desirable to compute the 
temperature and salinity contributions to the geostrophic flow sepa- 
rately. Furthermore, salt is not as easily exchanged between the air 
and sea as is heat, and therefore, it may be more acceptable to use 
historical data for salinity information than for temperature informa- 
tion. By directly computing the thermal and haline components of the 
geostrophic flow the significance of air-sea heat exchanges on thermo- 
haline flow may be more easily examined. 

To summarize, the purpose of this thesis is to study a simpli- 
fied method of computing geostrophic surface currents over large 
ocean areas that gives essential agreement with currents computed 
by the standard geostrophic method. The derived method must be 
computationally faster on digital computers than the standard method. 
Furthermore, that the method ought to express, at least to first 
order, the thermal and haline components separately. 
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JL THE GEOSTROP H1C CURRENT 



The Equations of Motion for Geostrophic Flow 



Fluid is said to be in geostrophic balance if the flow is unac- 
celerated and the horizontal pressure gradient is balanced by the 
Coriolis force. The simplified equations of motion for geostrophic 
flow in rectangular coordinates (x axis east, y axis north and z axis 
positive down from the sea surface) are: 



fkxV 



H 



Q V 



0) 



_ 3P 

g a 3Z 



( 2 ) 



w h ere 



f = 

oj - 



4> = 
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a 



v H p ^ 



2cj sin dp is the Coriolis parameter (l/ sec) 
the angular velocity of the earth (radians/sec) 
the geographic latitude (degrees) 
the horizontal velocity (cm/sec) 

2 

the local acceleration of gravity (cm/sec ) 
specific volume (cm /gm) 

2 

the horizontal pressure gradient (dynes/cm ) 



3P = 
3Z 



/ £ 

the vertical pressure gradient (dynes/cm ) 



A 

k ~ the vertical unit vector 
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Given the horizontal pressure gradient along a level surface the 
geostrophic flow along the surface could be determined directly from 
Equation 1. However, direct measurement of the horizontal pressure 
gradient in the sea has not been achieved. Several schemes have been 
developed to compute the horizontal pressure gradient from the mass 
distribution as determined from the measured temperature and 
salinity distributions. 

The Helland- Hansen Equation 

Geostrophic computations, also often referred to as dynamic 
computations or the dynamic method, are discussed in any elementary 
physical oceanography text, for example Neumann and Pierson ( 1967). 
Fomin (1964) discusses in considerable detail the application and the 
limitation of the dynamic method. Yao (1967) gives an excellent re- 
view of various computational schemes for computing geostrophic 
currents. What might be called the standard or classical method is 
based on an equation originally derived from Bjerknes* circulation 
theorem (1900) by Sandstrom and Helland- Hans en (1903), (as cited by 
Proudman, 1963). However, the resulting equation is commonly re- 
ferred to in the literature as the Helland- Hansen equation. The 
Helland - Hansen equation, written in terms of the horizontal gradient 
of the geopotential between two levels P^ and P^ is, 
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x ( v i ■ v 2 } _ v h § a ( x > y> p) c1 p 

p. 



( 3 ) 



w h e x* e 




x, y, p)dp - the geopotential of level P with respect 

/ 2 

to (cm/sec) 



P 



(Vj " V^) - the horizontal velocity at P^ relative to the 

level P^ (cm/sec) 



The geopotential difference between two levels in the ocean is 
determined by the integral in the right-hand term of Equation 3. 

The integral is numerically evaluated (using the trapazoidal rule) 
from computed vertical density distributions, determined from the 
measured vertical temperature and salinity structure. Using two 
stations the measured temperature and salinity structure can be used 
to compute the horizontal gradient in geopotential along an isobaric 
surface, and Equation 3 can be evaluated to give the geostrophic flow 
between the two stations. 



Use of Helland- Hansens 1 formula suffers from several limita- 
tions (Fomin, 1964). The feasibility of computing ocean currents 



The Limitations of the Dynamic Method 



using the dynamic method is based on the following assumptions: 
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1) The horizontal velocity and the horizontal pressure gradient 
is balanced by the Coriolis lorce. 

2 ) The horizontal velocity and the horizontal pressure gradient 
become negligible at a moderate depth below the sea surface. 

3) The field accelerations and frictional forces can be neglected 
(Svei drup, 1 947). 

That these assumptions are realistic for the large scale mo- 
tions, at least sometimes, has been illustrated by Stommel (1965). 
Even when the assumptions of geostrophy are satisfied, the neces- 
sary measurements are difficult to obtain to the desired accuracy. 

Errors in the Dynamic Method 

Errors in the dynamic method have been estimated by many 
authors, however, no completely satisfying error analysis of geo- 
strophic procedures is available. The errors in computing the geo- 
strophic current from standard hydrographic measurements (tem- 
perature and salinity at discrete levels) can be divided into the two 
categories, measurement errors and computational errors. 

Errors in the geostrophic current due to errors in the hydro- 
graphic measurements have been discussed by Wooster and Taft 
(1958), Reid (1959) Fomin (1964), and others. Following Reid ( 1959) 
the measurement errors in classical hydrographic work are sum- 



marized below: 
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1. Incorrect temperatures due to errors in temperature 
readings. 

2. Incorrect temperatures due to errors in the location of 
the measurement. 

3. Incorrect salinities due to errors in the salinity readings. 

4. Incorrect salinities due to errors in the location of the 
measurement. 

5. Incorrect temperature due to time and space variability. 

6. Incorrect salinity due to time and space variability. 

The errors in computing geostrophic currents due to computa- 
tional procedures have been discussed by Rattray (1961), Fomin 
(1964), Yao (1967) and others. These errors are summarized below; 

1. Errors in interpolation of point hydrographic data between 
sampling levels. 

2. Rounding and truncation errors in the numerical integra- 
tion. 

Due to the dependence of the measurement and computational 
errors on one another, no completely rigorous analysis is yet avail- 
able for the accumulated error in determining geostrophic currents 
by the dynamic method. The analysis by Fomin (1964) appears to be 
the most complete to date, but his discussion is not fully documented. 
The magnitude of the total error in dynamic height due to measure- 



ment and computational errors has been estimated independently by 
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several authors. Some of the available estimates are presented in 
Table 1. Error in the geostrophic current due to the indicated error 
in dynamic height is given in Table 1 for stations separated by 
100 km at 45° latitude. 



Table 1. Estimates of the error in the dynamic height and geo- 
strophic surface current relative to 1, 000 db. 



Author 


Dynamic Height 
Error (dynamic meters) 


V elocity 
Error (cm/sec) 


Wooster and Taft 


0 . 01 1 


1. 1 


(1958) 


Reed and Laird 


0. 003 


0. 3 


(1966) 


Kollmeyer (1964) 


0. 018 


l. 8 


Fomin ( 1964) 


0. 004 


0. 4 



The significance of these errors on geostrophic computations 
depends on the magnitude of the current. Certainly if the currents 
at mid-latitudes are only a few cm/sec then the relative error may 
be 50 to 100 percent. Fomin (1964) states in his conclusions con- 
cerning the accuracy of the geostrophic current computations that, 
"computational errors may completely distort the result. " Reid 
(I960) estimates that the combined error in salinity, temperature, 
pressure and position are on the order of 20 percent. While these 
relative errors seem prohibitive, an important point must be made 
here: from a practical point of view it is significant to know that 



1 1 

the geostrophic currents are weak. Therefore, if the currents are 
less than a few centimeters per second, this fact can be determined 
even if the relative error is large. 

Verification of Geostrophic Surface Currents 

A few direct current measurements have been collected in sup- 
port of geostrophic computations from hydrographic measurements. 
The classic case is a comparison reported by Wust ( 1924 ) of 
Pillsbury's current observations ( 1885 - 1889 ) in the Straits of Florida 
to geostrophic currents computed from hydrographic data collected 
in 1888 and 1912 . Observed currents are in remarkable agreement 
with the computed geostrophic currents, von Arx ( 1962 ) reports a 
comparison of currents observed with the geomagnetic electro- 
kinetograph (GEK) and dynamic topography showing substantial agree- 
ment in the Gulf Stream. Broida ( 1966 ) made simultaneous measure- 
ments of currents and hydrographic casts in the Florida Straits. 

His measured surface currents were in fair agreement with the com- 
puted geostrophic surface currents. Reed ( 1965 ) compared geo- 
strophic currents and currents measured by parachute drogue in the 
Alaska Stream and found excellent agreement between the two values 
for the surface current. Smith ( 1931 ) showed a definite association 



between the geopotential topography in the Grand Banks Region and 
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the drift paths of icebergs that established the basis for the U. S. 
Coast Guard hydrographic surveys of this region which continue 
today. 



IIL STANDARD GEOSTROPHIC COMPUTATIONS 



Prior to exploring attempts to simplify the standard geo- 
strophic computations these computations should be reviewed. 
Several schemes have been suggested for computing geostrophic 
currents: (1) the geopotential scheme using the Helland- Hansen 

formula (Equation 3), (2) the acceleration potential scheme of 
Montgomery and Stroup ( 1962), and (3) the isanosteric contour slope 
scheme of Werenskjold (1935, 1937). These three methods have 
been reviewed by Yao (1967). Only the geopotential scheme will be 
reviewed here, because of its role in later sections of this thesis, 
and its more frequent use over the other methods. 

The standard computations in the geopotential scheme as ex- 
pressed by the Helland- Hansen equation are expressed in the follow- 
ing equation: 



(V 



1 



-v 2 ) = 



< aD A - AD B> 

fAn 



(4) 



where 



AD a - AD 
A B 



AD 



= the difference in the anomaly of dynamic height 

/ 2 

at stations A and B (cm/sec) 

= Z6.AP., where 6\ is the average specific 
iii i 

volume anomaly in the i pressure interval 
(10 ^ cm *V gm) 
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(V 1 - V^) = the magnitude of the relative current normal 

to a line joining the two stations (cm/sec) 



In order to carry out geostrophic computations with Equation 4 
several complicated formulae must be used, or the empirical ex- 
pression of these formulae in tables (LaFond, 1951). The formulae 
necessary to calculate in situ specific volume anomaly from meas- 
ured or interpolated temperature and salinity are summarized by 
Yao (1967). 

Salinity as a function of chlorinity is given by Knudsen’s (1901) 
equation: 



It the measurement of electrical conductivity is to be used to deter- 
mine salinity or chlorinity then other formulae or tables must be 
used. 

Given the chlorinity the quantity cr can be computed. 



S = 0. 0 30 + 1. 8 05 Cl ( Knuds en, 1901) 



(5) 



where 



S - salinity (%o) 



Cl r: chlorinity (%o) 



o 



0- = -0. 069 + 0. 4708C1 - 0. 001570C1 2 + 0.0000398C1 3 

O 



( 6 ) 
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<r t = - T + (<r o + °* 1 324) [ ' * A T + B T (<J o ’ °- 1 324)1 (7) 



w here 



"t : ( Ps, T 0»‘° 3 



(Bjerknes and Sandstrom, 1910) 



and 



P 



S, T, o 



- the in situ density at atmospheric pressure 



o = a if T = 0° C 
o t 



r _ | (T - 3. 98) w , r T + 28 3 , 

~T L 503. 570 J L T + 67. 26 1 

A t = T(4. 78 67 - 0. 098 I 5T + 0. 0010843T 2 )10 

B t = (18. 0 30 - 0. 8 164T + 0. 0 1 667 T 2 ) I o" 6 
T = temperature (°C) 



The specific volume of a sea water sample in situ can then 
be computed from the equation of state (Ekman, 1908) 



a S, T, P a S, T, o " Pa S, T, O 



10* 9 { 



4886 



1 + 0. 0001 8 3P 



-[227 + 28. 33T (8) 



0. 551T 2 + 0. 004T 3 ] + P10 _4 [l055 + 9- 50T - 0. 158T 2 ] 



_ o cr - 28 _ 

1. 5P T 1 0 [ 147. 3 - 2. 72T + 0. 04T 



10 



_ a -28 - 

2.-, . . O ,4 r 



-P 1 0 _4 ( 32. 4 - 0. 87 T + 0. 002T )] +( -2^ ) [4. 5 



+ 0. IT - P10 _4 (l. 8-0. 6 T)] } 
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where 



ti 0 ^ = in situ specific volume (cm /gm) 
o, 1 > P 



a = specific volume at atmospheric pressure (cm /gm) 

S, T > o 

P = pressure (decibars) 

Frequent oceanographic practice uses the specific volume 
anomaly rather than the specific volume. 



6 a S, T, P ' a 35, O, P 



(9) 



where 



6 - specific volume anomaly (cm /gm) 



35, O, P 



- specific volume at S - 35%o# T - 0°C, and P 
decibars (cm /gm) 



the latter quantity is computed from the following equation: 



Pa 



[- 



4886 



35, O, P 35,0,0 35, O, O l 1 + 0. 0000183P 



- 227 



+ 0.01055P - (a _ - 28)( 14.73-0.0000324P)] 10 

3 D t o > u 



-9 



( 10 ) 



where 



a 35, o, O " °* 972643 cm /g m 



^s.o.o 328 - 126 



Specific volume is often expressed in the following way: 
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S. T 



+ 6 



S, P 



+ 6 



S, T, P 



(Sverdrup, 1933, cited (II) by (11) 
LaFond, 1951) 



where 



A 



S, 



T 



A 



S, 



T 



6 



S, 



P 



6 T, P 



6 



S, 



T, P 



- 6 + 6 + 6 and is called the thermosteric anomaly 

O 1 Dj 1 

(cm'Vgm) (Montgomery and Wooster, 1954 ) 

= 0.0273596 - 1 0 “ 3 cr t /(l - lO' 3 ^) 

~ Q S, T, O ' a 35, O, O 

" 6 S, O, P ■ 6 35, O, O 

" 6 35, T, P " 6 35, O, O 

= 6 • A S, T • 6 S, P " 6 T, P (6 S, T, P iS an ° rder ° f mag - 
nitude smaller than any of the other terms and normally 

neglected in geostrophic calculations. ) 



The individual terms in Equation 11 are given in tables (Sverdrup 
et ah , 1942 and LaFond, 1951). 

These formulae illustrate the computational complexity of 
standard geostrophic computations. To carry out these computations 
synoptically, on a hemispheric grid such as that used by the U. S. 
Navy Fleet Numerical Weather Central is not desirable if it can be 
avoided. The question remains can geostrophic computations be 
simplified, without giving up necessary accuracy? 
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IV. PAST EFFORTS TO SIMPLIFY GEOSTROPHIC COMPUTATIONS 

Several attempts have been made to simplify geostrophic com- 
putations. The attempts can be divided into three categories: those 

which (1) neglect the pressure dependent terms in the specific volume 
anomaly, (2) assume a simple relationship between the density of sea 
water and the variables (temperature, salinity, and pressure), and 
(3) determine the correlation function between two of the variables in 
particular water masses. The first approach is the most direct and 
will be the most generally applicable, and is discussed in the follow- 
ing section. 



Thermosteric Anomaly 

Montgomery and Wooster (1954) showed that the "thermo steric 

anomaly, "A , can be substituted for the specific volume anomaly 

b, I 

in geostrophic computations under some conditions without signifi- 
cant loss of accuracy. One condition is that the computations be 
limited to the upper layers, above 1, 000 db where the pressure 

terms (6 and 6 ) are not large; thus the thermosteric anomaly 

b> ir I , ir 

is a good approximation to the specific volume anomaly. 

Performing geostrophic computations on 47 hydrographic sta- 
tions from the Atlantic and Pacific oceans, using both the complete 
specific volume anomaly and the thermosteric anomaly, Montgomery 
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and Wooster showed that, except for one station in the Atlantic 
Ocean, the pressure terms contribute at most five percent to the 
station to station difference in the anomaly of dynamic height. These 
authors stated (page 66): 

The conclusion is reached that for hydrostatic computa- 
tion limited to the upper 500 db or l, 000 db, especially 
in the Pacific Ocean, if extreme precision is not required 
and if significant convenience is gained, the pressure 
terms may be neglected; in other words the thermosteric 
anomaly Ag ^ may be used in place of the complete specific 
volume anomaly 6. 

While this is clearly a simplification of the standard geo- 
strophic computations, the use of the thermosteric anomaly requires 
determination of c r and cr ^ Since the sole purpose for determin- 

ing these quantities is to evaluate the thermosteric anomaly in geo- 
strophic computations, a scheme which would allow geostrophic 
computations without the determination of these preliminary 
quantities would be more efficient. An interesting point made by 
Montgomery and Wooster is that for all practical purposes the pres- 
sure terms are nearly linear from the surface to the 2, 000 db level 
or deeper. 



Temperature-Salinity Correlation Methods 



One of the earliest attempts to simplify dynamic computations 
is that of Stommel (1947) through the use of temperature and 
salinity correlation. Since specific volume is a function of 
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temperature, salinity and pressure, all three of these variables must 
be accurately known in order to compute this quantity. However, if 
a correlation exists between any two of the independent variables, 
say temperature and salinity, then the specific volume anomaly could 
be expressed in terms of only two quantities, since the correlation 
function determines the relationship between the other two quantities. 
Since temperature as a function of depth is an easier and more com- 
mon measurement, Stommel suggested geostrophic computations 
could be made from temperature structure measurements alone 
using the temper ature - salinity correlation to determine the salinity; 
from this the specific volume anomaly and then the dynamic height 
could be evaluated. Oceanographers have long recognized that tem- 
perature and salinity are correlated in certain water masses 
(Sverdrup et al. , 1942). That is, for each temperature there will 
be a small range of salinity in a given water mass. 

Assuming that a satisfactory temperature - salinity correlation 
exists in a given water mass, then using this temperature - salinity 
correlation a new set of tables can be constructed giving the specific 
volume anomaly as a function of only two terms: 

6 = [T] + [T, P] (12) 

The first term on the right-hand side, [ T ] , is in reality the thermo- 
steric anomaly, p , which is a function of the temperature only 
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because of the established temperatur e - salinity correlation. The 
second term, [T,P], includes S Q _ and 6 but is fully deter- 

mined from temperature and depth in the given water mass because 
of the established temperature-salinity correlation. It is clear that 
the accuracy of the tempe rature - salinity correlation method is de- 
pendent on the nature of the tempe rature - salinity correlation. For 
each temperature there will be a certain finite range of salinity 
(R g ) which will represent the uncertainty in specifying the salinity 
from the temperature - salinity correlation. The uncertainty in 
salinity decreases with the slope of the tempe rature - salinity curve 
on the temperature - salinity diagram. Obviously this is due to the 
fact that for small slopes a small error in the temperature leads to 
a large change in the implied salinity. The uncertainty in salinity 
will also be large whenever seasonal variation and mixing lead to a 
poor temperature-salinity correlation. Values of R can be deter- 
mined at each level from the scatter of temperature - salinity pairs 
around the mean. Associated with the value of R at each level for 



a given water mass is an uncertainty in specific volume R^. Sum- 
mation of the values of R over the water column during the compu- 

o 

tation of dynamic height gives a measure of the uncertainty in the 
calculated dynamic height introduced by the uncertainty in the tem- 
perature-salinity correlation. 



Stommel applied the temperature - salinity correlation technique 
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to stations of the International Ice Patrol (3275 to 3278) off the Grand 

Banks, ATLANTIS stations 1637 to 1642 across the Gulf Stream, and 

some stations (unidentified) from the Sargasso Sea. He found that 

the use of temperature - salinity correlation was unsuited off the 

Grand Banks as the temperature-salinity correlation is poor in this 

region. His results are not surprising since the Grand Banks region 

is a mixing region for water masses of the Labrador Current and 

Gulf Stream (Kollmeyer, 1966). In the Gulf Stream he found that the 

temperature-salinity correlation method can be applied, but only if 

the stations are grouped, as the temper ature - salinity correlation 

changes across this current. However, in the Sargasso Sea where 

temperature-salinity correlation is excellent the method seems well 

suited. Stommel concludes (page 91): 

As a result it appears that in certain restricted regions 
the temperature- salinity diagram may be used for rough 
dynamic computations. For more details survey work 
where great accuracy of the results is desired the method 
using the tempe rature - salinity diagram is clearly un- 
suitable. 

LaFond (1949) applied the temperature - salinity correlation 
scheme in the Marshall Island region in an attempt to determined 
the geostrophic flow in this region using bathythermogram s. He 
used existing hydrographic data to establish the temperature- 
salinity correlation function for this region. LaFond's conclusion 
(page 236) summarizes the results of his studies: 
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To use bathythermogram s in determining relative cur- 
rents, there are several prerequisites which must be 
met: (1) the bathythermogram must extend at least to 
900 ft. , (2) the temperature-salinity relation must be 
established, with consideration of seasonal and geo- 
graphic effects, and (3) the effects of internal waves must 
be largely eliminated. If these conditions are attainable 
the direction of relative currents can be established from 
bathythermogram s. The results of this test indicate that 
the speed of the current (0/305 db) can be within 25 per- 
cent of those obtained from Nansen bottle and reversing 
thermometer data (O/l, 000 db). 

Yausi (1955) used temperature-chlorinity correlation in the 
Kuroshio Current, and later (1957) in the adjacent seas of Japan to 
determine the dynamic height anomaly from measurement of tem- 
perature alone. Approximating cr ^ by the following expression: 

0- t = A + BT + CT 2 + DCL (13) 



and using the temperature-chlorinity relationship to express 
chlorinity as a function of temperature for the water mass he de- 
veloped an approximate equation for the dynamic height anomaly. 



AD = 0. 016 - 0. 000241033 



J 

I 

0 



1 , 000 



TdZ 



+ 0. 000010461 



i 

I 



1 , 000 



T dZ - 0. 00000033426 






1,000 



3 

T dZ 



(U) 



Yausi further simplified this expression by determining the linear 
correlation function between the first integral, and the second and 
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third integral. Using the derived correlation function, Equation 9 
can be written: 

1 , 000 

AD= (0.2637 +0.018 2) +(0.000 l 543±0. 00000 15) J TdZ (15) 

0 

Comparing the dynamic height anomalies computed using 
Equations 14 and 15 to those computed by the standard method Yausi 
found standard deviations of 0. 0484 dynamic meters and 0. 0772 
dynamic meters respectively. These are at least four and seven 
times the accepted error in the standard method of dynamic compu- 
tations (Wooster and Taft, 1958). Computing the surface currents 
through a section extending south off Shionomisaki along approxi- 
mately 135. 7°E, Yausi found good agreement between the geo- 
strophic surface currents found from Equation 14 and those computed 
by the standard dynamic computations. 

Yausi attempted to achieve further simplification of Equation 
14 by expressing the dynamic height at the sea surface relative to 
1, 000 db in terms of the dynamic height of the surface relative to 
shallower levels, i. e. 300, 400, 500, and 600 db. Again the linear 
correlation function was determined as the relationship between the 
quantities. The error in determining ^ from the thermal 

structure above 600 db was 0. 048 dynamic meters, above 500 db was 
0. 051 dynamic meters, above 400 db was 0. 081 dynamic meters, 
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and above 300 db was 0. 129 dynamic meters. He concluded that a 
reasonably accurate measure of the geostrophic surface currents 
could be obtained by using only the thermal structure above 500 db. 

Yausi ( 1957) attempted to apply these methods to the adjacent 
seas of Japan, east of Honshu. While the patterns of dynamic 
topography obtained by the standard dynamic method and using the 
simplified equations of the form of Equation 14 and 15 were similar, 
an uncertainty of 0. 08 dynamic meters was inherent in the simple 
method due to poor temperature - chlorinity correlation. This is not 
unexpected because of the nature of the temperatur e-chlor inity cor- 
relation curve in the adjacent seas and the strong seasonal variability 
in this region. 

Attempts were made to apply the method over other parts of 
the North Pacific using stations 100 to 150 in the seventh cruise of 
the CARNEGIE. The difference between the dynamic height computed 
by Equation 14 and the standard method was less than 0. 06 dynamic 
meters, except for stations 126 and 131 off the west coast of the 
United States where the difference between the two methods reaches 
0. 18 dynamic meters. 

One concludes that the temperature-salinity or temperature- 
chlorinity methods are applicable in specific regions. However, it 
is not suited to other regions where the tempe rature - salinity curves 
have small slopes or in regions where seasonal variations and mixing 
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lead to a large scatter of temperature - salinity pairs. Furthermore, 
tempe rature - salinity correlations would have to be established for 
each water mass, and the correlation functions altered at water mass 
boundaries. Since water mass boundaries are not fixed such corre- 
lation will be difficult to determine. 



U. S. Navy Schemes 



Two groups in the U. S. Navy make ocean surface current com- 
putations on a synoptic basis over large ocean regions. One group 
uses techniques developed at the Naval Oceanographic Office 
(NAVOCEANO) principally for application in the Northwest Atlantic 
Ocean (James, 1966). The other group, Fleet Numerical Weather 
Central (FNWC), computes surface currents over the Northern 
Hemisphere every 12 hours (Hubert, 1964). Both groups rely on 
synoptic thermal structure data transmitted from ships at sea. The 
density and accuracy of thermal structure data has been discussed 
by Wolff (1964). 

NAVOCEANO Scheme 



NAVOCEANO techniques were originally reported by Gibson 
(1962). The technique is based on a relationship between the hori- 
zontal surface temperature gradient and the surface current speed. 
Discussing a hydrographic section along the 50th meridian Gibson 
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summarizes (page 4) the basis for the scheme: 

These data and sections for other ocean areas (not shown) 
form the basis for the analytical approach described 
below. Symmetrical undulation of the isotherms indicates 
four major water masses. 

Upon crossing each mass the surface current changes 
direction in an orderly manner; that is, the circulation 
is cyclonic for cold waters and anticyclonic for warm 
waters. There is also general agreement between the 
magnitude for temperature gradients and current velocity. 

/N 

If V is the surface current velocity, k a vertical vector, 
positive outwards, and AT is grad T, the relation, 

V = k cross AT holds in principle. This relation, 
analogous to that which applies for straight air flow sug- 
gests that water bands can be treated as greatly elongated 
air masses. 

James (1966) discusses the application of Gibsons 1 suggestion 
to synoptic analysis of ocean surface currents. Plotting observed 
currents against observed horizontal surface temperature gradients 
in the Northwest Atlantic, a set of curves is derived giving the sur- 
face current as a function of horizontal surface temperature 
gradient. The curves are shown in Figure l for the Gulf Stream, 
Sargasso and Labrador water masses for summer and winter condi- 
tions. The surface current is obtained by determining the sea sur- 
face temperature gradient and reading current speed from the ap- 
propriate curve. Direction of flow is assumed to be parallel to the 



isotherms. 



James states (page 60): 



CURRENT SPEED inmi/hr) 



ZB 




TEMPERATURE GRADIENT (°F/l5 nmi) 



Figure 1 



NAVOCEANO relationship between surface current speed and 
horizontal surface temperature gradient in the Northwest 
Atlantic (James, 1966). 



29 



This system, aside from the ease of computation, has two 
advantages: (1) use of input based on synoptic temperature 

data is apt to be more reliable than the use of climatological 
means, and (2) the direction of the flow is fairly accurate. 

Interstate Electronics Corporation (1968) made an evaluation 

of the geostrophic prediction techniques used by NAVOCEANO. In 

this study 7, 000 hydrographic station pairs were selected from the 

archives of the National Oceanographic Data Center. Only those 

station pairs were selected that were consecutive and separated in 

time less than 24 hours. Geostrophic computations were carried out 

for the 0 db and 125 db levels relative to the 1, 000 db level. The 

geostrophic currents at 0 db and 125 db between station pairs were 

correlated with the horizontal temperature gradient at these levels. 

The 125 db computation should have overcome local effects due to 

heat exchange across the sea surface. The results of the study were 

negative showing no correlation between the horizontal temperature 

gradient at the surface or at 125 db and the geostrophic current at 

these levels. 

NAVOCEANO determines the wind drift component of the sur- 
face current independently and adds this to the geostrophic com- 
ponent to obtain the total surface current. Other forcing factors are 
not considered. Wind drift is determined from curves relating the 
surface drift to wind speed, duration and fetch (James, 1966). 
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FNWC Scheme 

Hubert (1964) presents the equation used by FNWC to compute 
surface geostrophic flow over the Northern Hemisphere oceans on a 
synoptic basis. 

(V - V ) = V 't (16) 

' 1 2 fT H V ’ 

where 

AZ - the depth of the level of no motion 
V^T = the horizontal gradient of T 

T = K. T r + K-T„-, where K, and K 0 are "tuning 
1 sfc 2 200 1 2 

constants " 

T ^ and T ^ the surface and 200 meter temperatures respec- 
tively. 

Hubert (1964) does not develop Equation 16, nor can this in- 
vestigator find the development for application to the ocean. How- 
ever, Equation 16 is of the same form as the thermal wind equation 
(Haltiner and Martin, 1957). Development of the thermal wind equa- 
tion requires that there be a linear relationship between the tempera- 
ture and density, that is, the thermal wind is parallel to the mean 
virtual isotherms with low temperature on the left. Such a relation- 
ship cannot exist in the sea since the density of sea water is a nonlinear 
function of temperature, salinity and pressure. Hubert points out 
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that the coefficients K, and K can be adjusted in areas where 

1 2 

salinity gradients are known to be significant. 

The use of only the surface and 200 meter temperature fields 
cannot be justified over much of the ocean, as these two fields are 
not necessarily representative of the fields at other levels. Clearly, 
more of the water structure than just the sea surface temperature and 
200 meter temperature is necessary to make meaningful geostrophic 
computations. 

The wind drift is computed from Wittings (1909) formula 




where 

W = the current velocity (cm/sec) 

V = the 24 hour mean geostrophic wind speed (M/sec) 

- n wind factor 1 ' 

The factor is adjusted to account for the mass transport as- 

sociated with waves and the change of velocity with depth. For the 
surface wind drift is taken to be 4. 8, and the surface wind drift 

current is assumed to be parallel to the geostrophic wind. 

Both the geostrophic current and the wind drift current compu- 
tations are carried out on a 63 x 63 linear grid system on a polar 
stereographic projection over the Northern Hemisphere (Figure 2). 
The grid point separation is given by the following expression: 
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Figure 2. FNWC polar stereographic projection and grid 
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X 



[ l + sin 60° ] 
[ 1 + sin <j>~ 1 



x 200 



(18) 



The projection is true at 60° N where the grid spacing is 200 nautical 
miles. 

The wind drift and geostrophic surface current components 
(i and j) are computed independently at each grid point, and these 
components are added to arrive at the combined wind and geostrophic 
u and v components. From these components the current magni- 
tude and current direction are determined. The isotachs are con- 
toured resulting in a total transport (Figure 3). 

In order to obtain a single chart containing both speed and 
direction the components are used to obtain the vorticity of the cur- 
rent flow and the stream function ^ is obtained by a relaxation 
solution of Poisson’s equation (Hubert and Laevastu, 1965). 



2 _ 3v 9u 

V ^ 9x 9y 



(19) 



However, this stream function analysis is only applicable to non- 
divergent irrotational flow. The stream function field corresponding 
to the current transport field shown in Figure 3 is shown in Figure 4. 

Direct verification of FNWC computed surface currents has 
not been possible. Indirect verification of surface currents has been 
made through the verification of sea surface temperature analysis 
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Figure 3. FNWC current transport 
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Figure 4. FNWC stream function analysis 
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I 



(Hubert and Laevastu, 1965). Sea surface temperature analyses are 
made twice a day at FNWC over the Northern Hemisphere (Wolff, 
1964). Sea surface temperature changes computed from the air-sea 
heat exchange are subtracted from the analysis. If the residual cor- 
relates with the temperature advection field determined from the 
surface current analysis then the currents are assumed to be correct. 
If the residual does not correlate with the temperature advection 
field then the currents can be "tuned" to achieve agreement. It 
should be noted that the air-sea heat exchange formulae used in this 
procedure are generally empirical and have not been subjected to 
widespread rigid verification. This indirect verification procedure 
may be questioned seriously when used for quantitative results. 

Discussion 

The methods for determining geostrophic flow described in 
this section represent attempts to simplify geostrophic computations. 
Each model had some success but all suffer from certain defi- 
ciencies. The temperature- salinity correlation schemes of Stommel, 
LaFond, and Yausi greatly reduce the number of computations re- 
quired to determine geostrophic currents. These methods could 
also reduce sharply the field measurements, allowing geostrophic 
surface currents to be determined from temperature measurements 
alone. However, temperature - salinity correlation is a regional 
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parameter, varying from water mass to water mass; and in indi- 
vidual water masses it varies in both time and space. Furthermore, 
in regions of intense mixing such as the Kuroshio -Oyashio confluence 
( Tully, 1964), temperature-salinity curves are so variable that cor- 
relation approaches are not applicable. 

Use of the thermosteric anomaly v Q instead of the specific 

o» 1 

volume anomaly reduces the number of computations, or number of 
tables, that need to be interpolated in geostrophic computations. 

The error introduced by this simplification is acceptable if the com- 
putations are limited to the upper 1, 000 db. However, even in this 
simplified procedure the quantities o* and a ^ must be determined; 
since these are not used beyond the geostrophic computations, their 
determination represents unnecessary expenditure of computation 
time. 



NAVOCEANO correlation curves between current and surface 
temperature gradient may provide adequate synoptic current informa- 
tion in some regions where strong surface temperature gradients 
are dynamically sustained, such as in the Gulf Stream. Over the 
world oceans this is not the case; the surface temperature gradients 
are small and often are determined by the local heat exchange and 
mixing processes. In fact no definite relationship was found be- 
tween the surface current and the horizontal temperature gradient 
in 7, 000 selected hydrographic station pairs in the Northwest 
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Atlantic where the NAVOCEANO curves (Figure i) were derived 
(Interstate Electronics Corporation, 1968). 

The application of the thermal wind equation to the ocean has 
no foundation. The simplifying assumption that the density of sea 
water is a linear function of temperature cannot be accepted 
( Fofonoff, 1962). The adjustment of the coefficients used by FNWC 
on the basis of sea surface temperature analysis verification must 
be questioned. Since the sea surface temperature field is developed 
using the surface current field to compute the advected heat the use 
of sea surface temperature in surface current verification is not an 
independent verification. 

Geostrophic surface currents are most valuable if determined 
over a short period of time, say over a seven day period (Lenczyk, 
1964). This will be possible only when synoptic fields of tempera-' 
ture and salinity are available. While such fields are not available 
at the present time, a more scientifically sound scheme of comput- 
ing geostrophic currents must be available when these fields do be- 
come available. Any scheme used must yield surface currents that 
are in agreement with observed surface currents or at least those 
computed by the standard geostrophic method. Every effort should 
be made to verify by direct measurements any scheme of indirectly 
computing surface currents. 
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V. TEMPERATURE. SALINITY GRADIENT SCHEME 



Development of the T-S Gradient Equation 



Consider the Helland- Hansen equation, Equation 3, for geo- 
strophic flow where the horizontal gradient is expressed in differ- 
ential notation and n is perpendicular to (v^ - v^) in the horizontal 
plane. 




(n, P) dp] 



P 



1 



(3) 



Equation 3 can be rewritten: 

P 



(V 



i " V 2 ) = 7 \ 2 [ d / 

1 2 f J a(n ’ 



P)1 p dp 



( 20 ) 



where the subscript p indicates the bracket quantity is 
evaluated at constant pressure. 

Assuming that sea water can be regarded effectively as a binary 
fluid system whose specific volume is a function of the three inde- 
pendent variables, temperature, salinity and pressure. 

CL = f (T, s, P) (21) 

Then Equation 22 can be written in the following form carrying out 
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the indicated differential operation (Reid, 1959). 



( v _ v ) = i r 2 [ { ^) + ( i£) d P 

(V 1 * 2 f J 3T g ( P dn 9S T , P P 



(22) 



Note that there is no term representing the compressibility of sea 
water because the operation in brackets is carried out at a constant 
pressure. 

Equation 22 gives the geostrophic velocity in terms of the hori- 
zontal gradients of temperature and salinity; a quantity (8a/ 9T) 

b) x 

specifying the dependence of specific volume on temperature at a 

given salinity and pressure, and a quantity ( 9a/ 9S) specifying 

1 > r 

the dependence of specific volume on salinity at a given temperature 
and pressure. These latter quantitites become the coefficient of 
thermal compressibility and saline contraction if each is divided by 
the specific volume. Note that the compressibility of sea water 
enters Equation 22 indirectly through the dependence of these quanti- 
ties on pressure. 

The question then arises, as to whether or not simple expres- 
sions be found for (9a/9T) and (9a/9S) such that Equation 

Jr T , a 

22 represents a substantial simplification over standard geostrophic 

computations without significant loss in accuracy. Two possible ways 

of expressing ( 9a/9T) and (9a/9S) are: (1) to use 

b, Jr 1 9 Jr 

empirical data giving the specific volume as a function of tempera- 
ture, salinity and pressure, such as Newton and Kennedy (1965) or 
(2) to use one of the available equations of state, such as that of 
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Ekman (1908), that have been numerically fitted to the available 
empirical data. 



Determination of (3a/9T) and (9a/9S) 

O, P T , P 



Empirical Data 

The most widely used P-V-T (pressure, specific volume, 

temperature) data for sea water are based primarily on Ekman's 

(19 08 ) compression determinations. He measured the specific 

volume of a sample of sea water taken from 3, 000 m at a station off 

Portugal, at two different salinities; 31. 13 %oand38.83 %o , obtained 

by dilution and evaporation of the sample, and at three pressures; 

200, 400 and 600 bars. V-T data often used is that of Forsch, et al. 

(1902) for different salinities at atmospheric pressure. Forsch, 

et al. used a total of 24 samples collected entirely from the surface, 

mainly from the Baltic, North Sea and the North Atlantic Ocean. 

A more recent and extensive set of measurements are those 

of Newton and Kennedy ( 1965). They carried out measurements of 

specific volume for three salinities (31. 52, 34. 99 and 41. 03 %o) at 

temperatures from 0 to 25°C in 5° C steps, and at pressures from 

1 to 1, 000 bars in 100-bar steps. The precision of the measure- 

5 

ments is reported to be better than seven parts in 10 . Because of 
the time lapse of about six decades between the measurements of 
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Ekman and Forsch, and those of Kennedy and Newton, the latter 

measurements should reflect any advance in technique and apparatus 

in that interval. Furthermore, P-V-T-S data available prior to 

Newton and Kennedy is not sufficient to determine the quantities 

(3a/3T) and (3a/SS) over the range of temperature, salinity, 

b, P 1 > P 

and pressure of interest in the ocean without extensive interpolation 
of the data. However, Ekman's compressibility data is internally 
consistant to a remarkable degree (Eckart, 1958; Li, 1967). 

The dependence of specific volume on temperature, a(T) 

b, lr 

at fixed salinity and pressure, and the dependence of specific volume 
on salinity, a(S) , at fixed temperature and pressure are il- 

1 , ir 

lustrated in Figures 5 and 6 respectively (Newton and Kennedy, 1965). 

Figure 5 shows that a( T) is a nonlinear function over the range 

O, P 

of variables shown, and of practical interest in this study. However, 

a(T) c ^ is a continuously increasing function over these ranges, 
b, Jr 

Figure 6 shows that a(S) is nearly a linear function over the 

t , Jr 

range of variables of interest in this study. 

The quantities of interest, (3a/3T) and ( 3a/3S) , can 

S, P T, P 

be determined by direct differentiation of a(T) and a (S) 

T , P T, P 

if suitable expressions can be found for these functions. Polynomials 
of progressively higher degree (first through fifth) were fit in the 
least square sense to Newton and Kennedy’s data to determine 
a (T) s , p over the three salinities, at pressures of 1, 100, and 200 



SPECIFIC VOLUME (ml/gm) 
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Figure 5. The dependence of specific volume as a function 
of temperature at fixed salinity and pressure 
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1 bar 

1 bar 
1 bar 

100 bars 
ZOO bars 

ZOO bars 

ZOO bars 



Figure 6. The dependence of specific volume as a function 
of salinity at fixed temperature and pressure 
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bars. The fitting was performed on the Naval Postgraduate School 
IBM/ 360 d igitai computer using the program LSQPOL (Jordan and 
Vogel, 1961). This program computes the coefficients of t lie poly- 
nomial, an estimate of the error in the coefficients, and the standard 
deviation of the computed points from the fitted points. The standard 
deviation of the computed points was less than the precision of the 
original data points for a second degree polynomial. Therefore, 
a(T) can be expressed to the accuracy of the original data points 

O, X 

by an expression of the form: 



a( T) = a (S, P) + A(S, P)T + B(S, P)T (23) 

U) X O 



The values of a (S, P), A(S, P), and B(S, P) are given in Table 2. 

A similar procedure was followed to determine a(S) over the 

1 > P 



range of temperature and the pressure previously given for a(T) 



'S, P 



a(S) was expressed to the precision of the original data points 

1 > ir 

by a polynomial of first degree. 



a(S) T p - a o (T,P) - C(T, P)S (24) 

The values a (T, P) and C(T, P) are presented in Table 3. 

Examining the values in Tables 2 and 3 inconsistencies are 
noted for which no physical reason is available. For example, in 
Table 2 the values of A(30. 52 %o, 100 bars) and B(30. 52 %o» 1 00 bars) 
are less than the values of these coefficients at 100 bars and 
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Table 2. a Q (S, P), A(S, P) and B(S, P) as function of salinity and 

pressure in cgs units. (From Newton and Kennedy, 1965). 



Salinity 

(%o) 


Pressure 

(bars) 


a (S, P) 
o 


A(S, P) 

x 10 6 


B(S, P) 
x 10 7 


30. 52 


i 


0.9761 


48. 71 


50. 00 


30. 52 


100 


0. 9715 


89. 14 


40. 00 


30. 52 


200 


0. 9672 


99. 93 


40. 71 


34. 99 


1 


0. 9726 


67. 28 


47. 14 


34. 99 


100 


0. 968 2 


88. 21 


43. 57 


34. 99 


200 


0. 9639 


109. 91 


40. 00 


40. 03 


1 


0. 9680 


90. 78 


40. 71 


40. 03 


100 


0. 9637 


1 1 1. 71 


37. 14 


40. 03 


200 


0. 9595 


1 26. 07 


36. 43 



Table 3. a (T, P) and C(T, P) as a function of temperature and 



O' 

pressure 


in cgs units (From 


Newton and Kennedy, 


1965). 


Temperature 

(°C) 


Pressure 

(bars) 


a 0 (T, P) 


C(T, P) 
x 10 5 


0 


i 


0. 9996 


-77. 02 


0 


100 


0. 9942 


-74. 23 


0 


200 


0. 9895 


-73. 24 


5 


1 


0. 9993 


-75. 13 


5 


100 


0. 9943 


-73. 24 


5 


200 


0. 9895 


-71. 35 


10 


1 


0. 9997 


-74. 14 


10 


100 


0. 9948 


-7 2. 25 


10 


200 


0. 9904 


-71. 35 


15 


1 


1. 0000 


-7 2. 34 


15 


100 


0. 9955 


-71. 35 


15 


200 


0. 9910 


-70. 36 


20 


1 


1. 0014 


-73. 24 


20 


100 


0. 9965 


-70. 80 


20 


200 


0. 9920 


-69. 46 


25 


1 


1. 0025 


-72. 34 


25 


100 


0. 9977 


-70. 45 


25 


200 


0. 9934 


-69. 46 
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34. 99 % This is opposite to the trend in the other data. Similar 
inconsistencies are found in Table 3 for the C(S, P) coefficients. 
These inconsistencies are not unexpected. They are due to the fact 
that few data points are used and small errors in the data points fit 
lead to large errors in the coefficients. However, such incon- 
sistencies in the coefficients makes interpolation of coefficients 
from Newton and Kennedy's data to other temperatures, salinities 
and pressures impossible. The use of additional data points would 
overcome this difficulty. Additional data points smoothed to com- 
pensate experimental error can be computed from the equation of 
state for sea water. 

Equations of State for Sea Water 

Several equations of state have been suggested for sea water. 
The earliest equation (Equation 8) is that suggested by Ekman (1908) 
(Bjerknes and Sandstrom, 1910) for existing P-V-T-S data. 



a - a Q (l + up) 



(25) 



where 



a 



- specific volume, (ml/gm) 



a 



o 



- specific volume at atmospheric pressure 



p - pressure, in decibars 
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H = 10' 9 {[4886/(1 + 1. 83 x 10‘ 5 p)] - (227 + 28. 33T 

- 0. 551T 2 + 0. 004T 3 ) + 10“ 4 p (105. 5 + 9. 50T 

- 0. 158T 2 ) - 1. 5 x 10‘ 8 TP 2 - 1 0 ~ 1 (cr ^ 

- 28[(147. 3 - 2. 72T + 0. 04T 2 ) - 10 _4 p(32. 4 

- 0. 87 T +0. 02T 2 )] + 10" 2 (a- - 28) 2 [4. 5 - 0. 1 T 

o 

- 10" 4 p(l. 8-0. 06T)] } 

Ekman’s work would indicate that Equation 25 would be applicable 
over the following ranges in the variables: 

Temperature -2* C to 26° C 

Salinity , 31. 13 %o to 38. 5 3 %o 

Pressure 0 bars to 600 bars 

LaFond (1951), however, gives the range of application as: 
Temperature -2° C to 30* C 

Salinity 21 %o to 38 %o 

Pressure 0 bars to 1, 000 bars 

Eckart (1958) carefully studied the available P-V-T-S data 
for pure water and sea water. He concluded that the equation of 

- L 

state is represented to the accuracy of the available data (2 x 10 
ml/gm) by the Tumlirz equation. 
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(P + P )( a - q q ) = \ (26) 

where 

P = total pressure in atmospheres 

P = 5890 + 38 T - 0. 375T 2 + 3S 
o 

a ~ 0. 6980 
o 

\ = 1779. 5 + 1 1. 25T - 0. 0745T 2 - (3. 80 + 0. 01T)S 

Eckart indicates that Equation 26 is a satisfactory fit of the available 
data over the following range: 

Temperature 0° C to 40° C 

Salinity 0 %o to 40 %o 

Pressure 0 bars to 1, 000 bars 

Fofonoff (1962) compared the Ekman expression to the Tumlirz 
equation and found that the maximum disagreement between the two 
equations was less than 2x10^ ml/gm over nearly the entire 
range of salinity, temperature, and pressure in the sea. Only at 
unusually high ocean temperatures (greater than 29° C) with 
salinities of 36 %o did the disagreement reach 3x10^ ml/gm. 
However, while specific volumes computed by the two equations 
agree to within the accuracy of measurements, quantities derived 
from the equation of state, such as the coefficient of thermal expan- 
sion, are in serious disagreement. 

Li (1967) reviewed the available P-V-T-S data and suggested 
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the Tait-Gibson equation as an equation of state. The Tait-Gibson 
equation is given by Li as; 

a = a S) T , i * (1 ‘ S x 10 " 3 ) C L °g [frrf] < 27 ) 

where 

C - 0. 315 a _ . 

o, T, 1 

2 

B' 1 ' = (2670. 8 + 6. 890 56S) + ( 19. 39 - 0. 070 31 78S)T - 0. 223T 

Equation 27 is a satisfactory fit to existing P-V-T-S data over 
the following range of variables: 

Temperature 0° C to 20° C 

Salinity 30 %o to 40 %o 

Pressure (absolute) 1 bar to 100 bars 

Equation 27 gives results that are in agreement with measurements 
to the experimental error in the P-V-T-S data. The difference be- 
tween the Tait-Gibson equation and Ekman's equation for sea water 
of 35 %o t 0° C from 1 to 1, 000 bars is no more than 1x10 ml/gm. 
At atmospheric pressure the agreement between the density of sea 
water from Knudsen's tables (1901), in common usage in oceanog- 
raphy, and Equation 27 is less than 3x10^ gm/ml over the 
chlorinity range of 1 5 to 22 %o and temperature range of 0 to 20° C. 
Li concludes (page 2073): "Ekman's very involved equation of state 

of sea water is equivalent to the much simpler expression given 
here. 11 However, as the Tumlirz equation proposed by Eckart did 
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not give the same values as the Ekman equation for derived properties 
such as the coefficient of thermal expansion, the Tait-Gison equation 
gives again different values. 

The dependence of specific volume on temperature a(T) c 0 

O j x 

is a function of temperature, salinity and pressure according to the 

three proposed equations of state discussed above. In each case the 

functional relationships are markedly different. Furthermore, 

quantities derived from these expressions such as the coefficient of 

thermal expansion , l/a( a a /9T) , may differ significantly. 

o, F 

Fofonoff (1962) compared the coefficient of thermal expansion of sea 
water 35 %o salinity at atmospheric pressure computed from the 
Ekman and Eckart relationships. The results are given in Table 4. 



Table 4. A comparison of the coefficient of thermal expansion of 
sea water l/a(3a/3T)g p x 10^ at 35 %o salinity and 
atmospheric pressure computed from the Ekman and 



Eckart equations 


of state (from Fofonoff, 


1962). 


Temperature (° C) 


Ekman 


Eckart 


0 


52 


80 


5 


1 14 


121 


10 


167 


161 


15 


214 


201 


20 


256 


237 


25 


297 


274 


30 


335 


311 



While the coefficient of thermal expansion computed by the two 
equations are not significantly different at temperatures around 10° C 
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these differences increase significantly at higher and lower tempera- 
tures. 

Equation 22 requires the value of (3a/3T) _ before it can be 

used for geostrophic computations. To illustrate how this quantity 

changes between the three equations of state, Ekman, Eckart, and 

Tait-Gibson, the value of ( 3a/3T) as a function of pressure for 

o, F 

a salinity of 34 %o for temperatures 0, 10, 20, and 30° C is plotted 
in Figure 7. The three proposed equations of state yield values for 
(3a/8T) in closest agreement at temperatures near 10°C at pres- 

O, F 

sures near one atmosphere. However, the values diverge toward 
higher and lower temperatures and higher pressures. There is no 
clear cut way to specify which equation of state would yield the best 
value of (3a/3T) c . Furthermore, direct differentiation of the 

Ekman, Eckart or the Tait-Gibson equations leads to rather compli- 
cated expressions for (3a/3T) as well as (3a/3S) . Such 

b, F T , F 

complicated expressions would give Equation 22 no advantage over 
the Helland- Hansen equation, Equation 4. 

Fofonoff (1962) has examined the coefficients of saline con- 
traction, (3a/3S) , as computed from the Ekman and Eckart rela- 

1 t F 

tionships and found that they differ by less than one percent. There- 
fore, in geostrophic computations using Equation 22 there is little 

advantage in either expression for ( 3a/9S) . Thus, while the 

1 > F 

value of (3a/3T)g p differs significantly between equations of state 



1/ gm ° C) 
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I I " 1 I I I 1~~T— 1 I I 

0 20 40 60 80 100 120 140 1 60 180 200 

PRESSURE (bars) 

Figure 7. Variation of (0a /9 T) g p as a function of pressure for a 
salinity of 34% 0 accord’ing to the Ekman ( — - — )» 

Eckart ( ), Tait-Gibson ( ) equations of state 

and Equation 23 with coefficients derived from the 
Ekman equation of state ( ) 
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the value of ( 9a/ 3 S) „ does not. 

-L t -t 

Another approach to the determination of (3a/3T)g p and 

( 3a/ 3S) would be to compute specific volume at sufficient values 

T , P 

of temperature, salinity, and pressure using one of the equations of 
state, and fitting the resulting values with polynomials as was at- 



tempted with the Newton-Kennedy P-V-S-T data. This would over- 
come the difficulty encountered in directly fitting the empirical 
P-V-S-T data caused by experimental errors in the few data points 
as the equations of state would smooth individual inconsistent points. 
The question is which equation of state to use to compute the specific 
volume? There is little evidence that any of the expressions gives 
more reliable values of specific volume over the range of variables 
of interest in computing geostrophic currents in the upper few 
thousand meters of the ocean (Wilson and Bradley, 1968; Li, 1967). 

The Ekman equation (Equation 25) has been the principle equa- 
tion of state used by oceanographers (LaFond, 1951), Tables by 
Bjerknes and Sandstrom (1910), Sverdrup (1942) and LaFond (1951) 
are based on this expression. Since these tables have been and con- 
tinue to be so widely used in oceanography the Ekman equation of 

state is selected to determine a(T) and a(S) . 

S, P T , P 

Specific volume was computed using Equation 25 over the range 
of temperature, salinity, and pressure of interest in this study: 
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Temperature 0°Cto30° C 

Salinity 30 %o to 40 %o 

Pressure 0 bars to 200 bars 

Values of specific volume were computed for all possible combina- 
tions of the variables at intervals of 2° C for temperature, 2 % 0 for 
salinity and 10 bars for pressure. 

Second degree polynomials were fit to the specific volume as 
a function of temperature at fixed salinity and pressure to yield 

a(T) , and first degree polynomials were fit to specific volume as 
S, P 

a function of salinity at fixed temperature and pressure to give 
a(S) . Again the fitting was accomplished on the Naval Post- 

1 y x 

graduate School IBM/360 digital computer using the program LSQPOL. 

The values of the coefficients A(S, P) and B(S, P) as in Equation 23 

and C(T, P) as in Equation 24, are given in Appendix I as a function 

of their respective independent variables. The errors determined 

for these coefficients and the standard deviation of the computed 

versus the original data points is also given in Appendix I. 

The standard deviation of the goodness of fit of the computed 

to original data points, is always less than 2x10 ml/gm. Since 

Eckart (1958) contends specific volume is not known any better than 

2x10^ ml/gm the use of higher degree polynomials is not justified. 

Direct differ entiation of a(T) (Equation 23) with respect to tem- 

o, P 

perature yields the following expression: 
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(— ) = A(S, P) + 2B(S, P)T (28) 

9T S, P 

The values ( 3a/3T) from this equation were compared to the 

O, x 

values given by the Ekman, Eckart, and Tait-Gibson equations of 
state in Figure 7. The quadratic expression of a(T)g p leads to 
as good agreement with the commonly accepted Ekman equation of 
state as do either of the other proposed equations of state (Equation 
26 or Equation 27). 

In the same way differentiation of a(S) (Equation 24) with 

t » P 

respect to salinity yields: 



(— ) 
' 3S ; 



T, P 



C(T, P) 



(29) 



Linear interpolation of the computed values of the three coefficients 
leads to Figures 8, 9 and 10 for A(S, P), 2B(S, P) and C(T, P) 
respectively. 



Temperature - Salinity Gradient Scheme 



Recall Equation 22. 



(V 



V 2> 




ds, 

— dp 

dn P 



( 22 ) 



Substitution from Equation 28 and Equation 29 for (3a/9T) and 

I 

(3a/3S) T p respectively in Equation 22 yields the following 



PRESSURE (bars) 
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SALINITY (%c) 

Figure 8. A(S, P) as a function of salinity and pressure 
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simplified expression for computing geostrophic currents: 




+ 2B(S, P)T - C(T,P)^] dp (30) 



D 



Using Equation 30 geostrophic flow through a dynamic section 
can be computed without computing specific volume as required in 
the classical Helland- Hansen Equation, Only the measured tempera- 
ture and salinity gradients along isobaric surfaces, and the coeffi- 
cients A(S, P), 2B(S, P) and C(T, P) w'hich can be obtained from 
Appendix I or Figures 8, 9 and 10 respectively are required to com- 
pute the geostrophic flow through a section. Because of the clear 
relationship between the measured temperature and salinity gradients, 
and for lack of better, Equation 30 will be referred to as the Tem- 
perature-Salinity Gradient method or simply the T-S Gradient 
me thod. 

Computational Form of Temperature - Salinity Gradient Scheme 

Equation 30 is rewritten in terms of finite differences for 
computation: 




( 31 ) 
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where 

A. = the average value of A(S, P) over the i^ 1 pressure 
interval (cgs units) 

2B. r the average value of 2B(S, P) over the i 1 * 1 pressure 
interval (cgs units) 

C. = the average value of C(T, P) over the i^ 1 interval 
(cgs units) 

Al\ “ the average horizontal temperature difference in the 
i l pressure interval (°C) 

~ the average temperature in the i t ^ 1 pressure interval (°C) 

AS^ = the average horizontal salinity difference in the i** 1 
pressure interval (%o) 

2 

AP. = the i^h pressure interval (dynes/cm ) 

An = the horizontal separation of the stations (cm) 

Advantages of the Temperature- Salinity Gradient 
Scheme of Computing Geostrophic Currents 

If Equation 30 is to be used over the standard geostrophic 
scheme (Equation 4) it should represent some advantage. The ad- 
vantages realized in using Equation 30 are: (1) simplicity of compu- 

tation, (2) speed of computation, and (3) additional insight into the 
individual significance of the temperature and salinity structure in 
geostrophic flow in the ocean. Since the much simpler expressions 
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for a(T) and a(S) reproduce values of specific volume to 

the accuracy of the experimental measure of this quantity, these 
advantages would seem to be gained without loss of accuracy over 
the range of temperature, salinity and pressure considered in this 
study. Furthermore the method is completely general and can be 
used under any circumstances that the standard method can be ap- 
plied. 



The computational simplification of the T-S Gradient scheme is 
apparent since no density computations are needed. Neither the 
relatively involved tables commonly used to calculate specific volume 
(or specific volume anomaly), nor the equations of state previously 
discussed need to be used. Those properties that are directly meas- 
ured in the ocean (temperature and salinity as a function of depth) 
can be used to compute the geostrophic flow. Only the three coeffi- 
cients A(S, P), 2B(S, P) and C(T, P) are needed. 

The speed advantage of the T-S Gradient calculations on a 
digital computer was established on an IBM 360 digital computer. 

The time to compute the geostrophic flow through a section consist- 
ing of two stations, by the standard method and by the T-S Gradient 
method was measured. Only the actual time to carry out the mathe- 
matical computations was determined. The average time to com- 
pute the geostrophic flow for a single station pair was 0. 25 psec for 
the T-S Gradient method, and 1. 80 psec for the classical method. 
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Therefore, the T-S Gradient method is approximately seven times 
faster than the standard method using this digital computer. 

Another significant advantage of the T-S Gradient scheme is 
the potential for separating the thermal and haline contributions to 
the geostrophic flow. The interaction that takes place between the 
thermal and haline contributions to the flow is through A(S, P), 
B(S,P) and C(T,P). Therefore, at least to first order, the first two 
terms in Equation 30 can be considered the thermal contribution to 
the geostrophic flow, and the last term the haline contribution. 
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VI. APPLICATIONS OF THE T-S GRADIENT METHOD 

Introduce tion 

It has been shown that the T-S Gradient scheme is simpler to 
use than the standard geopotential scheme. It remains to evaluate 
the T-S Gradient scheme for making geostrophic computations over 
large ocean areas, say over the North Pacific. 

Clearly, additional advantage is gained in these calculations 
if the average value of the coefficients A(S, P), 2B(S, P) and C(T, P) 
relative to a selected reference level did not vary greatly over the 
ocean such that average values could be used and Equation 31 written: 



(V 



V.) = -d- [K.2T.AP. 
2 fAn 1 i l l 



+ K-ST.AT.AP - K.SaS.AP.] 



(32) 



where 



K , K and K - weighted mean values of A(S, P), B(S, P) and 

i j 

C(T,P) respectively between the surface and 
the selected reference level (cgs units). 

If the T-S Gradient scheme can be applied using Equation 32 then 
only the three fixed value coefficients and the measured temperature 
and salinity structure would be needed to make geostrophic computa- 
tions. To determine if Equation 32 is applicable, three questions 



must be answered: 
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1. What reference level should be used? 

2. Relative to the selected reference level, what is the 
spatial variation in the weighted mean values of the coeffi- 
cients (A(S, P), 2B(S,P) and C(T,P)? 

3. Relative to the selected reference level, what is the time 
variation of the weighted mean values of the coefficients 
A(S, P), 2B(S, P) and C(T, P)? 

Selection of the Reference Level 

The question of the proper reference level to select is open to 
question (Defant, 1961, and Fomin, 1964). Hopefully the level se- 
lected would be a level of little or no motion. Many reference 
levels have been used to compute surface currents (Neumann and 
Pierson, 1966). FNWC uses 200 m over the entire Northern Hemi- 
sphere oceans (Hubert and Laevastu, 1967). Reid (1958) used a 
reference level of 500 db in the California Current region, and 
l, 000 db reference level in his investigation (1961) of the dynamic 
topography in the entire Pacific Ocean. Reed and Laird (1966) used 
1, 000 db for dynamic computations in a section off central Cali- 
fornia. Reed and Taylor ( 1965) used 1, 000 db for geostrophic calcu- 
lations in the Alaska Stream. Stommel ( 1965) summarizes the 
problem of selecting a reference level (page 20): 
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Thus the choice of the reference level for geostrophic 
calculations becomes mostly a matter of taste, and we 
should admit that that is ultimately intolerable. The 
determination of the level of no motion is not a matter 
for debate, but for direct measurement. 

There are simply not sufficient direct current measurements in the 
oceans to select suitable reference levels. However, Fomin (1964) 
shows that it may be better to select a shallow reference level of 
low velocity than a substantially deeper level where the velocity may 
approach zero, thus reducing the computational errors. He points 
out that in some cases the selection of a deep reference level may 
completely distort the computed currents. 

For this study a level of 1, 000 db is selected as this is pos- 
sibly near an average value for the oceans. Furthermore, more 
hydrographic data is available over most of the ocean to 1 , 000 db 
than to deeper levels. Lastly, to go to a deeper reference level 
could introduce computational errors that may overcome any better 
approximation to no motion. Stommel (1965) states (page 19) 
‘'Fortunately, the choice of reference level has less effect upon 
velocity in the very surface layer than on those of deep water. M 



Determination of the Weighted Mean Value of K[, 
and K 3 Relative to 1, 000 db in the Summer and 
Winter for the Pacific Ocean 



Determination of 



K 2 and 



K. 



relative to any reference 



level requires that the distribution of temperature and salinity be 
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known to this level. These distributions are available from his- 
torical data deposited at oceanographic data centers. However, it 
requires a great quantity of data to accomplish oceanic coverage. 

As a first approximation, atlas values of temperature and salinity 
over the Pacific Ocean were used to determine weighted mean values 
of the coefficients relative to the 1, 000 db level. 

One of the most comprehensive atlases of temperature and 
salinity distributions of the entire Pacific Ocean is that by 
Muromotsev (1963). This atlas gives the fields of temperature and 
salinity at selected levels between the sea surface and the deepest 
part of the ocean. Fields are given for summer and winter condi- 
tions at the following levels: 0, 10, 25, 50, 100, 150, 200, 300, 

400 and 500 m. Below 500 m the seasonal variations are considered 
negligible and only single fields of temperature and salinity are 
available at each level: 600, 800, 1, 000, 1, 200, 1, 500 m, and 

thereafter every 500 m to the bottom. In all 30, 000 deep water 
hydrographic stations were used to compile the atlas supplemented 
in the upper layers by other temperature data. 

The values of K ^ and were determined for summer 

and winter conditions by interpolating the Muromotsev atlas tempera- 
ture and salinity fields at the following levels: 0, 10, 25, 50, 100, 

150, 200, 300, 400, 500, 600, 800 and 1, 000 m at the intersection 
of the lines of latitude and longitude given in Figure 11. The 



resulting weighted mean values of K^, and for summer and 

winter are shown in Figures 11 and 12, 13 and 1 4, and 1 5 and 16 



respectively. It can be seen from these Figures that the variation 
in the coefficients over the entire ocean, or from season to season 
does not change by more than a few percent. Using fixed values of 
Kj, and over the entire ocean relative to the selected refer- 

ence level of i, 000 db should lead to only a small percent error which 
is certainly acceptable in the geostrophic computation. The values of 
the coefficients selected are: 



Kj - 81. 0 x 10 ^ (cgs units) 

— -7 

“ 86. 5x10 (cgs units) 

— - 5 

- -74. 0x10 (cgs units) 

At most an error of five percent will be introduced into the geo- 
strophic surface current relative to 1, 000 db in the Pacific by the 
selection of these values. 



The Significance of Water Structure 

The use of fixed value coefficients relative to the 1, 000 db 
reference level in the Pacific appears to be justified. To further 
support this hypothesis, station pairs were selected in various re- 
gions of the world oceans and geostrophic calculations using the 
standard method and the T-S Gradient method using fixed values of 
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Figure 11. Variation of over the Pacific Ocean in the summer 
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Figure 12. Variation of over the Pacific Ocean in the winter 
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Figure 13. Variation of over the Pacific Ocean in the summer 
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Figure 14. Variation of over the Pacific Ocean in the winter 
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140 120 100 80 




L-1 



74.0 74.0 — 72.9 



73.0; — 73.0 — 73.2 — 73.3 — 73.5 — 73.6 




73.3 — 73 3 — 1 73 2 * — 73.2 73. 4 — 73. 5 — 73.5 



73.0 i; — 73.0 — 73.0 73.0 — 73.3 — 73.5 — 73.6 






73.8 — 73.9 — 73.9 — 73.9 — 74.0 




75.1 _ 75.1 75.1 — 75. 1 — 75.0 75.0 75.0 



120 140 160 180 160 140 120 100 80 
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60 



Figure 15. Variation of over the Pacific Ocean in the summer 
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Figure 16. Variation of over the Pacific Ocean in the winter 
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the coefficients. Lack of agreement would indicate that adjustment 
of the coefficients is required in some water masses. 

Stations used are shown in Table 5 where region, station 
identification, latitude and longitude, and the computed surface speed 
by each method are listed. The stations were selected to be approxi- 
mately perpendicular to the mean flow in each region. The complete 
velocity profile between the surface and the selected reference level 
shows small and generally not significant deviations between the two 
methods. The velocity profiles computed by both methods are given 
in Appendix II. The velocity profiles for each station pair computed 
by the two methods agree to better than ten percent at all depths. 

The comparison made in Table 5 shows the agreement required 
between geostrophic surface currents computed by the two methods. 
None of the values computed by the T-S Gradient method disagree by 
more than five percent from the values computed by the standard 
method. Therefore, the use of average coefficients in the T-S 
Gradient method to calculate geostrophic flow in the surface layer 
relative to the 1, 000 db reference level is justified. These results 
are encouraging. However, the need for direct verification of in- 
directly computed currents is great. Very few attempts to verify 
indirectly computed currents have been made, some are reviewed 
in Chapter II of this thesis. Clearly, more experiments are needed 
to establish the relationship between indirectly computed currents 



Table 5. Comparison of geostrophic surface currents computed by the standard method and the 
gradient method. 
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Table 5. Continued, 
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and actual observed currents. 

Verification of Indirectly Computed Currents 

Two field experiments were performed to investigate the ap- 
proximation of computed geostrophic surface currents to measured 
currents. These experiments are discussed in some detail in Ap- 
pendix V. The surface currents computed by the classical scheme, 
the T-S Gradient scheme and by FNWC were compared to the meas- 
ured surface currents. 

In both experiments the surface currents computed by the T-S 
Gradient scheme were identical with those computed by the classical 
scheme. In the first experiment the FNWC computed currents were 
in the opposite direction to the measured surface currents and the 
currents computed by both the T-S Gradient scheme and the standard 
geostrophic scheme. FNWC currents are heavily weighted toward 
the wind drift component. The computed geostrophic currents (4. 0 
to 9. 5 cm/sec) were only about half as great as the observed cur- 
rents (7. 1 to 14. 1 cm/sec) but they were in the correct direction. 

In the second experiment as in the first the computed currents 
and measured currents agree in direction, but the geostrophic com- 
puted currents were again slower (5. 4 cm/sec) than the measured 
currents (9. 7 cm/ sec to 16. 8 cm/sec). The FNWC computed cur- 
rents were the same speed ( 1 2. 7 cm/sec to 25. 5 cm/sec) and 
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direction as the measured currents. 

Thus in both experiments the geostrophic currents were nearly 
half the directly measured currents. This might be completely ex- 
plained by the fact that the mass structure was not sensed at the 
same time as the currents were measured. Even though these meas- 
urements were only separated by approximately 24 hours the mass 
structure may have changed or in fact never been in balance with the 
local currents. Furthermore, the geostrophic currents were calcu- 
lated from the average mass structure over at least a 48 hour period, 
and one would expect the average gradients to be less than the maxi- 
mum gradient. One fact these experiments show is the need for 
more field experiments to establish the validity of indirectly com- 
puted currents. 

In both experiments the T-S Gradient method and the standard 
geostrophic computations lead to essentially the same currents. 

That is, the simpler T-S Gradient computations give equivalent re- 
sults to the standard geostrophic computations and the advantages 
of simplicity are gained without giving up accuracy. Another ad- 
vantage of the T-S Gradient scheme is the separation of the thermal 
and saline contributions to the flow. An example of this advantage 



is presented in the next section. 
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The Relationship Between the Thermal Component of the Geostrophic 
Surface Velocity and the Total Geostrophic Surface Velocity 
in the Gulf Stream Water on the Grand Banks 

One of the advantages of the T-S Gradient method over standard 
geostrophic computations is the capability to separate the first order 
thermal and haline components of the geostrophic flow. To illustrate 
this feature of the method the thermal component and the total geo* 
strophic surface velocities were computed for 33 station pairs of 
the U. S. Coast Guard, International Ice Patrol stations in the Gulf 
Stream water (Kollmeyer, 1966) off the Grand Banks (Appendix IV). 
The linear correlation of the thermal component and the total geo- 
strophic velocity for the 33 station pairs was determined (Moynihan, 
1968 ): 

V = 0. 39V + 2. 84 
g t 

where 

V - the total geostrophic velocity (cm/sec) 
g 

V^ = the thermal component of the geostrophic velocity 
( cm/ sec) 

The standard deviation of the total geostrophic current from the 
linear correlation function between it and the thermal component is 
less than four (4) cm/sec. Therefore, in this water mass total geo- 
strophic currents might be determined using only temperature data 
through the relationship above. 
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VIL HEMISPHERIC COMPUTER APPLICATION 

Introduction 

The Fleet Numerical Weather Central, Monterey, California 
is the only facility attempting to determine surface currents on a 
hemispheric scale. The background of their computational technique 
has been given in an earlier section of this thesis (Chapter IV). For 
lack of any other hemispheric fields the FNWC 63 x 63 grid fields 
will be used in this study. The following fields are available through 
FNWC over their 63 x 63 grid: 

1. Synoptic temperature fields between the surface and 
1, 200 feet at 100 foot intervals. 

2. Climatological (monthly means) temperature fields be- 
tween the surface and 1, 200 feet at 100 foot intervals. 

3. Climatological (summer and winter) temperature and 
salinity fields in the Pacific and Atlantic, derived from 
published atlases (Muromtsev, 1963 and Bohnecke, 1936). 

Attempts were made in this study to use these fields to investigate 
the geostrophic flow through the T-S Gradient computations in the 
surface layers of the North Pacific. 
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Synoptic Fields 

During April, 1968, Equation 33 was programed on the Naval 
Postgraduate School IBM 360 digital computer to run over the FNWC 
synoptic temperature field for 0000 GMT, March 5, 1968. The pro- 

gram is found in Appendix III. While the program will calculate the 
salinity term in Equation 32, a constant salinity of 35 %o was as- 
sumed for lack of a salinity field. The output of the program is the 
i and j components of the velocity between the grid points, and the 
vector sum of these components at the southwest grid point of a four- 
point element as is shown in Figure 17. The direction of the total 
current was determined by hand computation from the ratio of the i, 
j components. Because of the large number of points (3889) in the 
63 x 63 grid only those in the North Pacific will be considered. How- 
ever, since the major oceanic circulation features are thought to be 
common to all oceans (Sverdrup, 1942) any conclusions drawn for 
the Pacific Ocean should be valid in the other oceans. 

For comparative purpose the FNWC Current Transport and 
Stream Function maps for the same period are shown in Figures 18 
and 19 respectively. Recall that these maps represent the sum of 
their geostrophic and wind driven components. Independent control 
for the geostrophic currents is provided by Tully (1964 after Reid, 
1961) who gives the surface geostrophic flow relative to 1, 000 db 
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Figure 17. Resolution of the i and j components of 
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(Figure 20). 

First, in Figure 19 the FNWC currents display many features 
of the established mean circulation in the North Pacific (Figure 20), 
i. e. , the subarctic gyre, subtropic gyre, and tropic gyre. Beginning 
at the Equator the equatorial system is shown, including the North 
Equatorial, South Equatorial, and Equatorial Counter Currents. 

The Kuroshio appears as a high speed flow and the Oyashio Current 
is present. The North Pacific Drift is apparent as is the separation 
of this flow into a northerly flow, providing the Gulf of Alaska circu- 
lation, and a southerly flow, the California Current. The magnitude 
of the recognized currents is reasonable, but in general the values 
are low in the narrow, concentrated streams such as the Kuroshio, 
and high in the regions of wide slow flows like the California Current. 

The currents computed with Equation 32 are plotted at selected 
grid points in Figures 18, 19, and 20. While some of the major cir- 
culation patterns can be identified, the current speeds are con- 
sistently low over the entire ocean. The computed currents are in 
poorest agreement in the western boundary area of the Kuroshio. 

Here Equation 33 over the 1, 200 foot temperature field gives a value 
which is low by an order of magnitude. The values in the North 
Pacific Drift are in agreement with those given by Tully (7 cm/sec 
to 10 cm/sec). The speeds in the Alaska stream also are low by an 
order of magnitude. The agreement is best in the California 
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Current. Some of the features of the equatorial current system ap- 
pear, but the speeds are high in the North Equatorial Current com- 
pared to those given by Tully. 

Several possible explanations are apparent for the low current 
speeds in the Kuroshio Current and in the Alaska Current. First 
these currents are narrow high speed flows. The grid spacing used 
in these calculations is so large as to completely straddle the cur- 
rents. Geostrophic computations assume that the velocity is uniform 
between any two stations in a section. The statement that any par- 
ticular stream is in geostrophic balance is meaningful only within 
the boundaries of the stream. If one station of a station pair is re- 
moved from the stream (and in many cases grid points used to calcu- 
late the velocity in a stream lie outside the stream boundaries) then 
geostrophic computations cannot be expected to be representative of 
the stream velocity. For this reason alone one would expect the geo- 
strophic velocity to be low in most cases since the mass gradients 
are rarely sensed at their maximum. Secondly, in both these flows 
a 1, 200 foot reference level is clearly too shallow. Lastly, only 
the thermal component of the current has been computed. Possibly 
the last two of these problems can be corrected by using the climato- 



logical temperature and salinity fields. 



89 



Climatological Fields 

To investigate the magnitude of the problems suggested in the 
previous section Equation 32 was programmed to run over the FNWC 
atlas climatological temperature and salinity fields on the FNWC 
CDC 3600 computer. The program is found in Appendix III. Since 
the atlas temperature fields were not digitized between the surface 
and 400 m, either the synoptic or the monthly climatological tem- 
perature fields must be used for temperatures above this depth. 

Thus fields from different sources are combined to make the calcu- 
lations. The compatability of the fields at any level and between 
levels has not been investigated. A reference level of 1, 000 m was 
selected as the average coefficients in Equation 32 had previously 
been determined for that reference level (Chapter VI). Furthermore, 
this was the reference level used by Reid (1961). The flow computed 
is shown in Figure 21 where isotachs indicate the surface current 
speed and arrows indicate the direction. 

The currents exhibit few of the features of the geostrophic cir- 
culation as given by Tully (1964), and many anomalies exist. In 
particular, there is a high transport region shown extending off the 
Japanese islands but the magnitude is too low by nearly half an order 
of magnitude. There is some indication of a westerly drift across 
the north Pacific, a southerly flow off the west coast of the United 
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Figure 21. T-S Gradient *tlas climatology 

surface geostrophic current field 
(cm/sec) relative to 1. 060 db 



91 



States and Mexico, and a cyclonic gyre in the Gulf of Alaska. The 
magnitude of the computed currents are in most reasonable agree- 
ment in the California Current region where values of a few centi- 
meters per second appear. In general the central portion of the 
subtropic gyre shows a strong northerly component which is not in 
agreement with known geostrophic circulation. Furthermore, the 
high transport in the western equatorial region is not in agreement 
with Tully. 



Remarks and Recommendations 



Remarks 

Equation 33 has been programed to run over the available 
63 x 63 fields of FNWC. The results of these computations have been 
compared to surface current computations presently made by FNWC 
and with geostrophic flow relative to 1, 000 m given by Tully (1964) 
for the Nortn Pacific. These comparisons show clearly that ac- 
ceptable geostrophic surface currents cannot be computed from the 
existing fields. While indications of accepted flow patterns appear, 
several limitations of the available fields are clear. 

For the synoptic fields one of the chief limitations is the fact 
that these fields only extend to 1, 200 feet. This is too shallow for 
a reference level over much of the ocean. Obviously, a second 
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limitation is the lack of salinity fields to complete the computations. 
Another limitation of the synoptic fields is the coarseness of the grid 
system. In many areas the grid spacing is larger than the known 
width of the existing currents. 

Using the existing climatological fields, the computations can 
be extended to 1, 000 db. Unfortunately, in the Pacific these fields 
do not represent adequately the mass field for geostrophic computa- 
tions. Derived from atlases, the fields are smoothed in an undeter- 
mined way. Furthermore, since the atlases present the temperature 
and salinity as horizontal fields at each level, the interpolation of 
these fields to the grid points undoubtedly will contribute errors. 

A further problem with the existing climatological fields is the re- 
quirement to use monthly climatological temperature fields in the 
upper 400 m along with atlas salinity fields, without establishing the 
compatability of the two fields. Again the coarseness of the available 
63 x 63 grid is not appropriate for geostrophic computations in the 
narrow current regions. If these computations are to be made then 
several adjustments must be made in the available fields. 

Clearly the available fields and coarse grid restrict geo- 
strophic computations over the ocean by any method. The more ac- 
ceptable agreement between currents computed by FNWC and 
"accepted 1 ' currents is due to their capability to "tune" the flow. 

This same procedure could be used to adjust the geostrophic flow 
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computed by the T-S Gradient scheme using the synoptic temperature 
fields which appear from this study to be the best available fields. 
However, "tuning*' the ocean to an arbitrary standard seems ques- 
tionable. Geostrophic computations should be objective and the use 
of arbitrary "tuning" factors removes this objectivity. Several 
recommendations are obvious which should be considered as possible 
improvements, hopefully leading eventually to the capability of com- 
puting actual geostrophic flow over the world oceans. 

Recommendations 

First, the 63 x 63 Northern Hemisphere grid used by FNWC 
gives a grid spacing of 37 3 nmi at the equator, 200 nmi at 60° N 
latitude, and 186 nmi at the pole. This grid spacing is too large for 
geostrophic computations over most of the oceans. The minimum 
spacing suggested by Reed and Laird (1965) for 160°W and 175°W 
longitude is an indication of more acceptable spacing. Their values 
are shown in Table 6. The basis for their recommendations is that 
the dynamic height difference should be five times the error in 
dynamic height computations. However, a grid spacing too much 
larger than this will lead to poor results. Since the dynamic height 
represents an integral of the water structure a grid should be de- 
veloped for geostrophic computations using the dynamic height 
gradient as one criterion in constructing a new grid and temperature 
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and salinity fields. 



Table 6. Minimum station spacing along 160° W and 
North Pacific (from Reed and Laird, 1965). 


175° W in the 




Optimum Station Spa 


cing (km) 


Current 


along approx. 
160° W 


along approx. 
175° W 


Alaska Stream 


16 


8 


Subarctic 


183 


110 


West Wind Drift 


61 


39 


Gyre 


> 137 


> 137 


Kuroshio Extension 


55 


34 



Another justification for altering the grid spacing considers 
the distribution of data. The smallest grid spacing is found at the 
poles where the lowest data density is available, while in the tropical 
regions the grid spacing is largest although this is where data is 
relatively more abundant (Wolff, 1964). Therefore, where the most 
data are available and the smallest scale features could be shown in 
the data, the grid scale is coarsest and these features are lost. 

New fields must be constructed if geostrophic computations 
are to be made. The best approach to the construction of these new 
fields would be to work with actual hydrographic data available from 
depositories of these data i. e. , the National Oceanographic Data 
Center, Washington, D. C* The requirements of geostrophic computa- 
tions should be used as a guideline in construction of new fields. 
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VIIL SUMMARY AND CONCLUSIONS 

This study has examined schemes for simplifying geostrophic 
computations. Several schemes have been previously proposed to 
simplify these computations, but except for the use of the thermo- 
steric anomaly proposed by Montgomery and Wooster (1954) none are 
satisfactory over the entire ocean. 

A simplification is proposed here which represents a signifi- 
cant simplification of dynamic computations without loss of accuracy. 
The scheme expresses the geostrophic velocity as three terms in- 
volving the observed variables, temperature and salinity. Only three 
coefficients are needed. These coefficients have been determined 
and are given in graphical and tabular form. If a fixed reference 
level of 1, 000 db is used it has been shown that fixed values of the 
coefficients can be used for geostrophic surface current computations 
over the entire Pacific Ocean without introducing more than a five 
percent error. The name T^S Gradient scheme is proposed for the 
simplification. It provides the following advantages over other 
schemes: 

1. Computations show no loss in accuracy from the classical 
scheme for most practical purposes. 

2. The computation of geostrophic currents is substantially 



simplified. 
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3. The scheme is general and can be used in any region of 
the ocean where the temperature, salinity and pressure fall within 
the range for which the coefficients have been determined. 

4. The scheme requires only the observed variables, tem- 
perature, salinity and pressure, and the coefficients A(S, P), 

2B(S, P) and C(T, P). 

5. The first order thermal and haline structure contributions 
to the geostrophic flow can be determined independently. 

6. The scheme is computationally faster than the classical 
scheme (seven times faster on the IBM 360). 

7. The scheme is particularly well suited to use onboard 
ships with small digital computers as involved equations of state 
are not used. 

The proposed scheme was used to compute surface geostrophic 
flow over the 63 x 63 FNWC northern hemisphere fields, the only 
hemispheric fields available at this time. The resulting computed 
surface currents were compared to surface currents in the North 
Pacific Ocean computed by FNWC, and by Tully (1964) from Reid’s 
(1961) investigation of the geopotential topography of the 0 db level 
relative to the 1000 db level. These comparisons lead to the con- 
clusions that the FNWC 63 x 63 fields are unsuited to geostrophic 
computations for the following reasons; 

1. The grid spacing is too large even at its smallest spacing. 
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2. The atlas derived fields do not adequately depict the mass 
distribution in the ocean. 

3. The synoptic temperature fields do not extend to deep 
enough level to reasonably approximate level of little or no motion. 
The computations over the 63 x 63 fields do show, however, that it 
is practical to make actual geostrophic computations over such a 
large ocean area on a synoptic basis if the desired fields are avail- 
able. 

Two proposals are made to improve geostrophic computations 
over hemispheric grids. First, that a new grid shouldbe constructed 
in which the grid spacing is compatible with geostrophic computa- 
tional requirements; second, new temperature and salinity fields be 
constructed from the available hydrographic data to replace the atlas 
climatological fields presently available. 

The T-S Gradient method can be used to compute independently 
the first order thermal and haline contributions to the geostrophic 
flow. It has been shown that the correlation of the first order 
thermal geostrophic surface current (V^) with the geostrophic sur- 
face current (V ) in the Gulf Stream water off the Grand Banks can 
g 

be expressed by a simple linear relationship. This relationship is 
useful for the determination of the total geostrophic surface current 
from temperature measurements alone. 

The reliability of the T-S Gradient scheme was tested by 
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comparing T-S Gradient and standard geostrophic computations in 
several differing water masses. The current velocity structure was 
nearly exactly the same in both cases (Appendix II). Never were the 
computations in disagreement by more than ten percent even though 
fixed coefficients were used. 

Two field tests were performed to investigate the approxima- 
tion of computed geostrophic surface currents (Appendix IV) to meas- 
ured surface currents. The currents were measured prior to the 
hydrographic measurements using parachute drogues. Geostrophic 
currents were computed from the measured temperature and salinity 
by both the classical geostrophic method and the T-S Gradient 
method, and FNWC currents were obtained from their Current 
Transport and Stream Function maps for the periods of the experi- 
ments. On the second experiment the geostrophic computations 
were also made with the reversing thermometer temperatures re- 
placed by the expendable bathythermograph temperatures. 

The geostrophic currents computed by both methods were 
identical for both experiments. The substitution of expendable 
bathythermograph temperatures raised the computed surface velocity 
by both methods approximately one cm/sec. The computed geo- 
strophic surface velocity was approximately half the observed cur- 
rent in both experiments. The FNWC currents were in the opposite 
direction to the observed flow in the first experiment but agreed 
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with the observed flow in the second. 

Certainly two experiments should not be considered as verifi- 
cation or condemnation of any scheme. The results of the experi- 
ments do show the need for additional experimentation before a great 
deal of confidence can be placed in any indirectly computed currents. 

It is proposed that this research be continued to further in- 
vestigate the relationships between indirectly computed surface cur- 
rents and observed currents (see Appendix V). 
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APPENDIX I 



Tables of A(S, P), B(S,P) and C(T,P) 
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6 6 
Table AI-1. a (S, P), 10 A(S, P), the error in A(S,P)[lO « ], 

°Y y 

10 B(S, P), the error in B(S, P )[ 1 0 « ] and the 

^ B(o, P ) 

standard deviation of 10 a(T) as a function of salinity 

o,P 

(%o) and pressure (db). (cgs units) 



S (%o) 


A(S. P) 


f A 


B(S, P) 


f B 


Sigma 


30 


56. 09 


P 

2. 06 


= 0 db 

46. 97 


0. 66 


2. 00 


31 


58. 89 


2. 04 


46. 55 


0. 66 


1. 99 


32 


61. 65 


2. 0 3 


46. 14 


0. 65 


1. 97 


33 


64. 41 


2. 01 


45. 72 


0. 65 


1. 95 


34 


67. 14 


2. 00 


45. 31 


0. 64 


1- 95 


35 


69. 85 


1. 99 


44. 90 


0. 64 


1. 94 


36 


72. 55 


1. 98 


44. 49 


0. 64 


1. 92 


37 


75. 23 


1. 97 


44. 08 


0. 63 


1. 92 


38 


77. 89 


1. 96 


43. 67 


0. 63 


1. 91 


39 


80. 53 


1. 96 


43. 27 


0. 63 


1. 90 


40 


8 3. 16 


1. 95 


42. 87 


0. 63 


1. 90 


30 


58. 72 


P = 
2. 03 


100 db 
46. 58 


0. 65 


1. 98 


31 


61. 49 


2. 01 


46. 16 


0. 65 


1. 96 


32 


64. 23 


2. 00 


45. 75 


0. 64 


1. 95 


33 


66. 96 


1. 99 


45. 34 


0. 64 


1. 93 


34 


69. 67 


1. 97 


44. 93 


0. 64 


1. 92 


35 


7 2. 36 


1. 96 


44. 52 


0. 63 


1. 91 


36 


75. 03 


1. 95 


44. 1 1 


0. 63 


1. 90 


37 


77. 69 


1. 94 


43. 71 


0. 62 


1. 89 


38 


80. 32 


1. 94 


43. 31 


0. 62 


1. 88 


39 


8 2. 94 


1. 93 


42. 90 


0. 62 


1. 88 


40 


8 5. 55 


1. 92 


42. 50 


0. 62 


1.87 



J% 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 
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Continued. 



A(S, P) 


€ a 


B(S, P) 


€ B 


Sigma 


61. 33 


P = 

2. 01 


200 db 
46. 19 


0. 64 


1. 95 


64. 07 


1. 99 


45. 77 


0. 64 


1. 93 


66. 80 


1. 97 


45. 37 


0. 63 


1. 92 


69. 50 


1. 96 


44. 96 


0. 63 


1. 91 


72. 18 


1. 95 


44. 55 


0. 63 


1. 89 


74. 8 5 


1. 94 


44. 15 


0. 62 


1. 88 


77. 50 


1. 93 


43. 74 


0. 62 


1.87 


8 0. 13 


1. 92 


43. 34 


0. 6 2 


1. 86 


8 2. 74 


1. 91 


42. 94 


0. 61 


1. 86 


8 5. 34 


1. 90 


42. 55 


0. 61 


1.85 


8 7. 91 


1. 90 


42. 15 


0. 61 


1. 85 


63. 93 


P = 

1. 98 


300 db 
45. 80 


0. 64 


1. 92 


66. 64 


1. 96 


45. 40 


0. 63 


1. 91 


69. 34 


1. 95 


44. 98 


0. 63 


1. 89 


7 2. 01 


1. 93 


44. 58 


0. 6 2 


1. 88 


74. 68 


1. 92 


44. 18 


0. 62 


1. 87 


77. 32 


1. 91 


43. 78 


0. 61 


1. 86 


79. 94 


1. 90 


43. 38 


0. 61 


1. 84 


8 2. 55 


1. 89 


42. 98 


0. 61 


1. 84 


8 5. 14 


1. 88 


42. 59 


0. 61 


1.83 


87. 71 


1. 88 


42. 19 


0. 61 


1. 82 


90. 27 


1. 87 


41. 80 


0. 60 


1.82 
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Table Al-1. Continued 



S (%o) 


A(S, P) 


f A 


B(S, P) 


f B 


Sigma 


30 


66. 50 


P = 
1. 95 


400 db 
45. 41 


0. 63 


1. 90 


31 


69. 19 


1. 94 


45. 01 


0. 6 2 


1. 88 


32 


71. 86 


1. 92 


44. 60 


0. 62 


1. 87 


33 


74. 52 


1. 90 


44. 21 


0. 61 


1. 85 


34 


77. 15 


1. 89 


43. 81 


0. 61 


1. 84 


35 


79. 77 


1. 88 


43. 41 


0. 60 


1. 83 


36 


8 2. 37 


1.87 


43. 01 


0. 60 


1.82 


37 


84. 95 


1. 86 


42. 62 


0. 60 


1.81 


38 


87. 52 


1. 86 


4 2. 23 


0. 60 


1. 80 


39 


90. 07 


1. 85 


41. 84 


0. 60 


1. 80 


40 


9 2. 60 


1. 85 


41. 45 


0. 59 


1. 79 


30 


69. 05 


P = 
1. 92 


500 db 
45. 03 


0. 6 2 


1. 87 


31 


71. 72 


1. 91 


44. 63 


0. 61 


1.85 


32 


74. 37 


1. 89 


44. 23 


0. 61 


1. 84 


33 


77. 00 


1. 88 


43. 84 


0. 60 


1. 82 


34 


79. 61 


1. 86 


43. 44 


0. 60 


1. 81 


35 


8 2. 20 


1.85 


43. 05 


0. 60 


1. 80 


36 


84. 78 


1. 84 


42. 65 


0. 60 


1. 79 


37 


8 7. 34 


1. 84 


42. 26 


0. 59 


1. 78 


38 


89. 88 


1. 83 


41. 88 


0. 59 


1. 78 


39 


9 2. 41 


1. 82 


41. 49 


0. 59 


1. 77 


40 


94. 9 2 


1. 82 


41. 10 


0. 59 


1. 77 
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Table AI-1. Continued. 



s (%o) 


A(S, P) 


( A 


B(S, P) 


€ B 


Sigma 


30 


71. 59 


P = 
1. 90 


600 db 
44. 65 


0. 61 


1. 84 


31 


74. 23 


1. 88 


44. 26 


0. 60 


1.83 


32 


76. 8 5 


1. 86 


43. 86 


0. 60 


1. 81 


33 


79. 46 


1. 85 


43. 47 


0. 60 


1. 80 


34 


8 2. 05 


1. 84 


43. 08 


0. 59 


1. 79 


35 


84. 62 


1. 83 


42. 69 


0. 59 


1. 78 


36 


87. 17 


1. 82 


42. 30 


0. 58 


1. 77 


37 


89. 71 


1.81 


41. 91 


0. 58 


1. 76 


38 


9 2. 22 


1. 80 


41. 53 


0. 58 


1. 75 


39 


94. 73 


1. 80 


41. 14 


0. 58 


1. 75 


40 


97. 22 


1. 79 


40. 76 


0. 58 


1. 74 


30 


74. 01 


P = 
1.87 


700 db 
44. 28 


0. 60 


1. 82 


31 


76. 72 


1.85 


43. 89 


0. 60 


1. 80 


32 


79. 32 


1. 84 


43. 49 


0. 59 


1. 79 


33 


81. 91 


1.82 


43. 10 


0. 59 


1. 77 


34 


8 4. 47 


1. 81 


42. 7 2 


0. 58 


1. 76 


35 


87. 02 


1. 80 


42. 32 


0. 58 


1. 75 


36 


89. 55 


1. 79 


41. 94 


0. 58 


1. 74 


37 


9 2. 06 


1. 78 


41. 56 


0. 57 


1. 73 


38 


94. 56 


1. 77 


41. 18 


0. 57 


1. 72 


39 


97. 04 


1. 77 


40. 80 


0. 57 


1. 72 


40 


99. 50 


1. 76 


40. 42 


0. 57 


1. 72 



Table AI-1. Continued. 



1 1 1 



S (Too) 


A(S, P) 


f A 


B(S, P) 


( B 


Sigma 


30 


76. 60 


P = 
1. 84 


800 db 
43. 91 


0. 59 


1. 79 


31 


79. 20 


1.82 


43. 52 


0. 59 


1. 77 


32 


81. 77 


1. 81 


43. 13 


0. 58 


1. 76 


33 


84. 33 


1. 80 


4 2. 74 


0. 57 


1. 75 


34 


86. 87 


1. 78 


4 2. 36 


0. 57 


1. 73 


35 


89. 40 


1. 77 


41. 94 


0. 57 


1. 72 


36 


91. 90 


1. 76 


41. 59 


0. 57 


1. 71 


37 


94. 40 


1. 75 


41. 21 


0. 56 


1. 70 


38 


96. 87 


1. 75 


40. 83 


0. 56 


1. 70 


39 


99. 32 


1. 74 


40. 45 


0. 56 


1. 69 


40 


101. 77 


1. 74 


40. 08 


0. 56 


1. 69 


30 


79. 08 


P = 

1. 81 


900 db 
4 3. 54 


0. 58 


1. 76 


31 


81. 65 


1. 80 


43. 15 


0. 58 


1. 75 


32 


84. 21 


1. 79 


42. 76 


0. 57 


1. 73 


33 


8 6. 74 


1. 77 


4 2. 38 


0. 57 


1. 72 


34 


89. 26 


1. 76 


42. 00 


0. 56 


1. 71 


35 


91. 76 


1. 74 


41. 62 


0. 56 


1. 70 


36 


94. 24 


1. 74 


41. 24 


0. 56 


1. 69 


37 


96. 71 


1. 73 


40. 87 


0. 56 


1. 68 


38 


99. 16 


1. 72 


40. 49 


0. 55 


1. 67 


39 


101. 60 


1. 72 


40. 1 2 


0. 55 


1. 67 


40 


104. 02 


1. 71 


39. 75 


0. 55 


1. 66 
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Table AI-1. Continued. 



S (%o) 


A(S, P) 


<A 


B(S, P) 


f B 


Sigma 






P = 


1000 db 






30 


81. 55 


1. 79 


43. 18 


0. 58 


1. 74 


31 


84. 09 


1. 77 


42. 79 


0. 57 


1. 72 


32 


86. 62 


1. 76 


4 2. 40 


0. 56 


1. 71 


33 


89. 13 


1. 74 


42. 03 


0. 56 


1. 69 


34 


91. 63 


1. 73 


41. 65 


0. 56 


1. 68 


35 


94. 01 


1. 72 


41. 27 


0. 55 


1. 67 


36 


96. 57 


1. 71 


40. 90 


0. 55 


1. 66 


37 


99. 01 


1. 70 


40. 53 


0. 55 


1. 65 


38 


101. 44 


1. 69 


40. 1 5 


0. 54 


1. 65 


39 


103. 85 


1. 69 


39. 78 


0. 54 


1. 64 


40 


106. 25 


1. 68 


39. 41 


0. 54 


1. 64 






P = 


1100 db 






30 


8 3. 99 


1. 76 


42. 81 


0. 57 


1. 71 


31 


86. 51 


1. 74 


42. 43 


0. 56 


1. 69 


32 


89. 02 


1. 73 


4 2. 06 


0. 56 


1. 68 


33 


91. 51 


1. 71 


41. 68 


0. 55 


1. 67 


34 


93. 98 


1. 70 


41. 30 


0. 55 


1. 65 


35 


96. 43 


1. 69 


40. 93 


0. 54 


1. 64 


36 


98. 87 


1. 69 


40. 56 


0. 54 


1. 64 


37 


101. 29 


1. 67 


40. 19 


0. 54 


1. 62 


38 


103. 70 


1. 67 


39. 8 2 


0. 54 


1. 62 


39 


106. 09 


1. 66 


39. 45 


0. 53 


1. 61 



108. 46 1. 66 39. 09 0. 53 1. 61 



40 



113 



Table AI-1. Continued. 



S (%o) 


A(S, P) 


f A 


B(S, P) 


€ B 


Sigma 


30 


86. 41 


P = 
1. 73 


1200 db 
4 2. 46 


0. 56 


1. 68 


31 


88. 91 


1. 72 


4 2. 08 


0. 55 


1. 67 


32 


91. 40 


1. 70 


41. 70 


0. 55 


1. 65 


33 


93. 86 


1. 69 


41. 33 


0. 54 


1. 64 


34 


96. 31 


1. 67 


40. 96 


0. 54 


1. 63 


35 


98. 74 


1. 66 


40. 59 


0. 53 


1. 61 


36 


101. 16 


1. 65 


40. 22 


0. 53 


1. 61 


37 


103. 56 


1. 65 


39. 85 


0. 53 


1. 60 


38 


105. 94 


1. 64 


39. 49 


0. 53 


1. 60 


39 


108. 31 


1. 63 


39. 12 


0. 53 


1. 59 


40 


1 10. 66 


1. 63 


38. 76 


0. 52 


1. 58 


30 


88. 82 


P = 
1. 70 


1300 db 
4 2. 10 


0. 55 


1. 66 


31 


91. 29 


1. 69 


41. 73 


0. 54 


1. 64 


32 


93. 76 


1. 67 


41. 35 


0. 54 


1. 63 


33 


96. 20 


1. 66 


40. 98 


0. 53 

• 


1. 61 


34 


98. 63 


1. 65 


40. 61 


0. 53 


1. 60 


35 


1 01. 04 


1. 64 


40. 25 


0. 53 


1. 60 


36 


103. 43 


1. 63 


39. 88 


0. 52 


1. 58 


37 


105. 81 


1.62 


39. 52 


0. 52 


1. 57 


38 


108. 17 


1. 61 


39. 16 


0. 52 


1. 57 


39 


110. 51 


1. 61 


38. 80 


0. 52 


1. 56 


40 


1 1 2. 84 


1. 60 


38. 44 


0. 52 


1. 56 
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Table 


AI-1. Continued. 










S (%o) 


A(S, P) 


( A 


B(S, P) 


€ B 


Sigma 


30 


91. 21 


P = 
1. 68 


1400 db 
41. 75 


0. 54 


1. 63 


31 


9 3. 66 


1. 66 


41. 38 


0. 53 


1. 61 


32 


96. 10 


1. 65 


41. 01 


0. 53 


1. 60 


33 


98. 52 


1. 63 


40. 64 


0. 52 


1. 59 


34 


100. 92 


1. 62 


40. 28 


0. 52 


1. 57 


35 


103. 31 


1. 61 


39. 91 


0. 52 


1. 56 


36 


105. 68 


1. 60 


39. 55 


0. 51 


1. 55 


37 


108. 04 


1. 59 


39. 19 


0. 51 


1. 55 


38 


110. 38 


1. 58 


38. 83 


0. 51 


1. 54 


39 


1 1 2. 70 


1. 58 


38. 47 


0. 51 


1. 54 


40 


115. 01 


1. 58 


38. 12 


0. 51 


1. 53 


30 


9 3. 58 


P = 
1. 65 


1500 db 
41. 40 


0. 53 


1. 60 


31 


96. 01 


1. 63 


41. 03 


0. 52 


1. 59 


32 


98. 43 


1. 62 


40. 67 


0. 52 


1. 57 


33 


100. 82. 


1. 60 


40. 30 


0. 52 


1. 56 


34 


103. 20 


1. 59 


39. 94 


0. 51 


1. 55 


35 


105. 57 


1. 58 


39. 58 


0. 51 


1. 54 


36 


107. 92 


1. 57 


39. 22 


0. 51 


1. 53 


37 


110. 25 


1. 56 


38. 86 


0. 50 


1. 52 


39 


112. 57 


1. 56 


38. 51 


0. 50 


1. 51 


39 


114. 8 7 


1. 55 


38. 15 


0. 50 


1. 51 


40 


117. 16 


1. 55 


37. 80 


0. 50 


1. 51 
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Table Al- 


1. Continued. 










S (%o) 


A(S, P) 


f A 


B(S, P) 


€ b 


Sigma 


30 


95. 93 


P = 
1. 62 


1600 db 
41. 06 


0. 52 


1. 58 


31 


98. 34 


1. 61 


40. 69 


0. 5 2 


1. 56 


32 


100. 73 


1. 59 


40. 33 


0. 51 


1. 55 


33 


103. 11 


1. 58 


39. 97 


0. 51 


1. 53 


34 


105. 47 


1. 57 


39. 61 


0. 50 


1. 52 


35 


107. 81 


1. 56 


39. 25 


0. 50 


1. 51 


3b 


110. 14 


1. 54 


38. 89 


0. 50 


1. 50 


37 


1 1 2. 49 


1. 54 


38. 54 


0. 49 


1. 49 


38 


114. 74 


1. 53 


38. 19 


0. 49 


1. 49 


39 


1 17. 03 


1. 53 


37. 84 


0. 49 


1. 48 


40 


119. 29 


1. 52 


37. 49 


0. 49 


1. 48 


30 


98. 27 


P = 
1. 60 


1700 db 
40. 71 


51. 32 


1. 55 


31 


100. 65 


1. 58 


40. 35 


50. 79 


1. 54 


32 


103. 02 


1. 56 


39. 99 


50. 30 


1. 52 


33 


105. 38 


1. 55 


39. 63 


49. 86 


1. 51 


34 


107. 71 


1. 54 


39. 28 


49. 47 


1. 50 


35 


1 10. 04 


1. 53 


38. 92 


49. 13 


1. 49 


36 


1 1 2. 34 


1. 52 


38. 57 


48. 8 2 


1. 48 


37 


1 14. 63 


1. 51 


38. 22 


48. 57 


1. 47 


38 


116. 90 


1. 50 


37. 87 


48. 36 


1. 46 


39 


119. 16 


1. 50 


37. 52 


48. 19 


1. 46 


40 


1 21. 41 


1. 49 


37. 17 


48. 07 


1. 45 
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Table Al- . 


L. Continued, 










S (%o) 


A(S, P) 


€ a 


B(S, P) 


€ b 


Sigma 






P = 


1800 db 






30 


100. 59 


1. 57 


40. 38 


0. 50 


1. 52 


31 


102. 95 


1. 55 


40. 02 


0. 50 


1. 51 


32 


105. 30 


1. 54 


39. 66 


0. 49 


1. 49 


33 


107. 63 


1. 52 


39. 31 


0. 49 


1. 48 


34 


109. 94 


1. 51 


38. 95 


0. 49 


1. 47 


35 


1 1 2. 24 


1. 50 


38. 60 


0. 48 


1. 46 


36 


1 14. 5 3 


1. 49 


38. 25 


0. 48 


1. 45 


37 


116. 79 


1. 48 


37. 90 


0. 48 


1. 44 


38 


119. 05 


1. 48 


37. 56 


0. 47 


1. 43 


39 


1 21. 29 


1. 47 


37. 21 


0. 47 


1. 43 


40 


123. 51 


1. 47 


36. 87 


0. 47 


1. 43 






P = 


1900 db 






30 


102. 89 


1. 54 


40. 04 


0. 50 


1. 50 


31 


105. 23 


1. 53 


39. 69 


0. 49 


1. 48 


32 


107. 55 


1. 51 


39. 33 


0. 49 


1. 47 


33 


109. 86 


1. 50 


38. 98 


0. 48 


1. 45 


34 


112. 16 


1. 48 


38. 63 


0. 48 


1. 44 


35 


1 14. 4 3 


1. 47 


38. 28 


0. 47 


1. 43 


36 


116. 69 


1. 46 


37. 93 


0. 47 


1. 42 


37 


1 18. 94 


1. 46 


37. 59 


0. 47 


1. 42 


38 


121. 17 


1. 45 


37. 25 


0. 47 


1. 41 


39 


1 23. 39 


1. 45 


36. 90 


0. 46 


1. 40 


40 


1 25. 59 


1. 44 


36. 56 


0. 46 


1. 40 
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Table AI-1. Continued. 



S (%o) A(S, P) € _A B(S, P) Sigma 

P = 2000 db 



30 


105. 17 


31 


107. 49 


32 


109. 79 


33 


1 1 2. 08 


34 


1 14. 35 


35 


116. 6 1 


36 


1 18. 85 


37 


1 21. 07 


38 


1 23. 28 


39 


125.48 


40 


127. 66 



1. 51 


39. 71 


1. 50 


39. 36 


1. 48 


39. 00 


1. 47 


38. 65 


1. 46 


38. 31 


1. 45 


37. 96 


1. 44 


37. 6 2 


1. 43 


37. 28 


1. 42 


36. 94 


1. 42 


36. 60 


1. 41 


36. 26 



0. 49 


1. 47 


0. 48 


1. 46 


0. 48 


1. 44 


0. 47 


1. 42 


0. 47 


1. 42 


0. 47 


1. 41 


0. 46 


1. 40 


0. 46 


1. 39 


0. 46 


1. 38 


0. 46 


1. 38 


0. 45 


1. 37 



Table AI-2. 10 C(T, P), 10“ c D v and the standard deviation of 10 a(S) T as a function of 

l 1 1 a ) J- » a 

temperature (°C) and pressure (db). 
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Table AI-2. Continued, 
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Table AI-2. Continued. 
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APPENDIX II 



Geostrophic Velocity Profile Comparison 
Standard Versus T-S Gradient Method for 
Various Regions in the Oceans 
(Stations Positions are Given in Table 5) 



Note: Original scales reduced in reproduction. 
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APPENDIX III 



Computer Programs for Hemispheric Computations 
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Program 1. A program to interpolate hydrographic data 
on the IBM/360 computer 



SUBRQLTINF LGINTP <N,V,D,SO,CV,NN) 

DI MFNSIGN D( 1( f >) , V< 1000) ,CV< ICAO ) , SCM 1-00 ) 
1 i 1 DC 186 J = i,50«. 

112 DO 186 I = 1 , N 

I F ( Si) ( J ) - D ( N ) ) 113,115,190 

115 CV( J ) = V( N ) 

Of) TO 190 

113 I F ( SD(J)-D(I )> 114,114,116 

114 CV< J ) =V( I ) 

GO TO 188 

116 IF< S D ( J * - D < 1 + 1>H2C ,118,186 
113 CV(J1=V(1 + H 

GO TO 188 

120 I F ( I - 1 ) 132,132,126 

126 XA = (SC<J>-D(m*(SD(J>-D<Uin*V<I-l)/ 

1( (DC I — 1 ) - D ( I ) ) * ( D ( I — 1 > — D ( I + 1 ) ) ) 

XB = ( SIX J»-D( 1-1) l*<SD( J)-D< 1 + 1) )*V< 11/ 

1 ( < D ( I )-()( I-lll»(D(ll-On+ll)l 
xc = ( SC ( J )-!.)( [-1II»(S0(JI-DIIII*VII*1)/ 
ii(nmiH)(Hii*(Dioii-n(im 
AN $IJ= X A + XB ♦ XC 

132 IF <( I + 2)-N) 133 , 133 , 134 

133 YA=< SD( vl » - 0 < I ♦ 1 > )*(SD<J>-0< I + 2) )*V( I )/ 

1 ( ( 0( I )-0( 1 + 1) ) *■ ( D ( I )-D( 1 + 2) ) ) 

YB = ( SD( J » -D ( 1 n*(SP(J)-D( 1 + 2) )*V( I ♦ 1 ) / 

u(n(i + ii-o(iii*ioii + n-o(i + 2m 

YC = (SD< ,))-()( I ) )*<SD(J)-D( I + l> i*V( 1 + 2)/ 

i(( 0 (I + 2 >-D(in*( 0 (I+ 2 )-D(I + im 

4NSD= YA+ YB+YC 

134 Zl--(St(JI-D(Hin*V(I)/ID(II-nil + ll] 

Z B = ( S C ( J ) -0 ( I ) > " V < 1+1) /(D( 1+1 )-D( I > ) 
ANSL=ZA+ZH 

I F ( 1-1)136,136,133 
136 CV( J) =< ANSO+ANSL ) /2. 

OLL = ( ANSD+ANSL + ANSi » / 3. 

GO TO 188 

138 IF< ( I+2I-N) 14C , 140,142 
140 C V ( J>=( ANSU+ANSD+ ANSI) /3. 

UD = ( ANSU+ANSD) /2. 

GO TO 188 

142 CV( J) =( ANSU + ANSD /2. 

ULL = ( ANSU + ANSL+ANSL ) /3. 

GO TO 188 
106 CONTINUE 
188 CONTINUE 
J = 5 0 1 
19C NN= J 
RETURN 
END 
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Program 2. A program to compute specific volume 
anomaly from hydrographic data on the 
IBM/360 computer 
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SUBROUTINE SGTSVA ( T * StDtSGT . SV.SVA ) 

ST— H < T-3.98>**2I/50 3.57I*((T*283. >/<T*67.26M 
CL = < S-. >30 1/1. fi )5 

SC = -» v_o9-4-1.4706*CL-.00157*CL*-*2*i.98E-5*CL**3 
4T=T*<4. 7867-.r98185*T*.OClC8‘t3*T**2)*l.E-3 
3 T = T * ( 18.03v>~. 8164*T*«C 1667*T ** 2 ) * 1 • F-6 
SGT=ST*( SO*. 1324»*< 1. -AT* BT* ( SO-. 1 324 > ) 

AF ST= 1 • /( l.+SGT*l.E-3l 
A=C*AFST*l.E-9 
8=4886. /< l.+1.83E-5*D> 

C=227.*28.33*T-. 55 1 * T **2 * . U04*T**3 
E =0* 1. E-4 
G = < SO-28. 1 /10. 

H=147. 3-2. 72*T+. 04* T* *2 
U=lo5.5+9. 5* T-. 158*T**2 
V= 1. 5*D**2*T*l.E-8 
W=32.A-.87“T*.( 2*T**2 
X = 4. 5-, 1 * T 
Y=l. 8-. J6*T 

SV=AFST-A*( B-C*E*<J-V-G*( H-E*W) *G**2*< X-E*Y ) ) 

4Z =• 972643 

YA = -227.*.0105f'*D 

YB=.C126*( 147. 3-. OG 324* D I 

AP=AZ-l)*AZ*(B+YA-YR|*l.E-9 

5 VA=SV-AP 

RETURN 

END 
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Program 3. A program to compute geostrophic cur- 
rents by the Helland- Hansen equation 
and the T-S Gradient equation on the 
IBM/360 computer 
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K L A l *8 IMSG 

REAL *4 I T l ( 1 2 ) 

DIMENSION SD < 5( * ,DA<5. 1 , TA(5t I ,SA< 50 ) , OKI 5'' I ,TB( 5-' I , SB ( 5«. » 

0 I V E N $ I ON TCA(£\.C) »SCA(50('( ,TCB(5lU ) , S C 4 ( 5' 0 1 , PL T ( 501 ) , DL S ( 5C0) 
DIMENSION HR T ( 500 ) .ANSI 500) . VELD ( 51.0 ) • A$ VA ( 5”t I , BSVAI 500 » 

JIMENSI ON JSVAI 500) .TSVAI 500 ) .CSVAI 50 j I .DOS ( 5L0 ) , VFLF ( 5( 0 ) 
DIMENSION IMSGI8) ,S0X<500) 

"0L I VALENCE ( IMSO( 1 I , I TL( L |» 

DATA ITK9I/' •/ , ITK IC I / •MIDDELftU*/, ITLI 1 1 )/*PG •/ 

OA TA I TL( 12 I / *P OX 85 ' / 

DATA 14/' •/ 

C .,.,0;76 .l.oCCCHS - . 0 L 0 7 3 6 

C 11 format ( 2E 5. 2 ,2 I 3 , 6 A8 I 
11 FORMAT (2F5. 1 ,2I3,8A3I 

14 FORMAT (8F1C.3) 

15 FORMAT (4F2C.C) 

00 43 1 =1 , 5C^ 

43 SDII I =1C0l'. 

READ (5,15) A , P , C 

44 R c AD (5,15) ALNM.STEP 
IF ( A l N M I 3L\,,3 T,45 

45 SO(l 1=3. 

DC 4'' 1=2, 5h j 
SOI I ) = SD ( 1-1 1 ♦STEP 
IF ( SD( I )-AlNM| 4 j,42 ,42 
40 CONTINUE 
42 N \ = I 

READ (5,111 DIS.PSI ,NA,NB,IMSG 
5o DO 57 1=1, NA 

READ (5,14) DAI I ) , TA( I ) ,SA( l ) 

IF ( T A ( I ) - 32 . )52 ,71 ,71 

71 T A < I (=( TA< I )-32. I *. 55555 
52 IF( SA ( I 1-25. I 55 ,55 ,57 

55 SA( I l=C..03+l. 8( 5*SA( I ) 

57 CONTINUE 

DO 55 I = 1 . NB 

READ (5,14) DB ( I ) , TB ( I ) , SB ( l ) 

IF ( TB< I )- 32. I 54,72,72 

72 TB( I )=( TB( I )-32. I*. 55555 

54 IF < SB ( I ) - 2 5. I £8,58,59 

58 S8( I ) = < • . C 3 ♦ 1 • 8( 5*Sti( I I 

55 CONTINUE 

CALL LGINTP ( NA ,TA ,DA ,SD ,TCA,NNA I 
CALL LGINTP (NA , S A , OA , SD , SC A , NN A ) 

CALL LGINTP (NP,T8,CB,SC,TCB,NNR) 

CALL LGINTP (NB,SH,CB,SC,SCP,NMi) 
i»2 = l. 45 8E-4 

PSJ=( 2. *3. 141 5^2 65/360. ) *PSI 
SPSI= SIMPSJ) 

IF ( SPSl.LT.D.l)SPSI=J.l 
X2= 1. /( 42*SPSI»DI S) 

DO 62 1=1, NN 

DL T ( I ( = TC A ( I ) -TCB ( I ) 

DL S ( I 1= SC A ( I ) - SCB ( I ) 

BR T ( I | = < TC A ( I I 4 TC FJ ( I ) l*C.5 
62 CONTINUE 
ND=NN- 1 
DO o4 1=1, ND 

SOL T = ( DL T ( I I +DL T ( I ♦ l ) I *C. 5 
DDLS =(DLS( I l+DLSI I ♦ 1 1 1*1.5 
BTDT = PDLT*( (BRT(1 ) ♦ BR T ( I ♦ 1 ) I*. 5) 

SD SC =SD( I +1 ) -SOI I ) 

ABDLT = A*BDLT* SDSD 
3BTDT = B*BTDT»$OSO 
CtiDLS = C*8DLS*SDSD 
SVEL = ( AB0LT4RHTDT4CBDLS ) *X2 
64 ANS (I)= U*BOl T*B*BTDT4C*BDLS)*(SD( I + D-SDI I ) ) 

VELDI NN ) =0, 

DO 66 1=2, NN 
JD = NN- ( 1-1 ) 

66 VELDI JD)=VtLD( JO+1 l+ANSI JO)*XZ 
DO 82 1=1, NN 
T S= TC A ( I ) 
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OS=SD( I I 

CALL SGTSV* ( TS ,SS .ns ,SGT ,SV, SVA) 

^ a s v a ( n =sva 
*<♦ i=i,nn 

T S= f Lt<( I ) 

SS*SCP( I ) 

;s = snt 1 1 

CUt SGTSV4 < TS,SS,DS,SGT,SV,SVA) 

^ \ 5VA< I ) = $ VA 

)i" 3 2 I =1 , NN 

?2 1SVM I )^= A 3 Vm ( I ) - 6 S V A ( I ) 

y' 9^ i =1 , no 

T 3 V A ( I)=(D5VA( l)>OSVA(Ul))*v.5 
4 C S V A ( I > - T S V A ( I ) * ( S D ( I ♦ 1 ) - $ D ( I ) ) 

30 $< MM) = l. 

S6 I=2,NN 
K =NlN- ( I -1 ) 

4 4 oos(k i=nos( k*i ^csva(k) 

DO 96 I = 1 t N N 

A3 VC LP( I ) = UOS< 1MX2 

WRITE <^>,9) IMjr, 

WRITE (6.1.1 

1J cormat ( 1 X . • D t f TH TEMPERATURE: SALINITY* , KX f ' A + 2R-C VELCCI 

1 T Y D DYNAMIC VELCCITY P*/* A 0 A B*, 

2 3 X f 'INTER GAL CM/StC HEIGHT CM/SEC 1 //) 

) ' i C 2 1=1. NO 

1.? *KfTE (6,10 SC ( I ) ,TCA( I ) ,TCB< I) ,SCA( I ) ,SCB< I ) ,VELO( I),DDS( I ) , 
l VC L P ( I ) .ANSI I ) 

U form at (1X.F5.L ,2F7. 2,2FR.3tl6X,lP3tl2.V4'X,lPFl?.3) 

*R I TL (6,13) SP( NN) ,TCA( NN) ,TC8( NN) ,SCA( NN) , SCR (NN 1 , VELD( NM ) , 

1 ;nS( inn ) , vflp ( Nf ) 

13 cr )R M A T ( lx f F5,i ,2F7.2,2F8.3,16X,1P^E1?.3/// ) 

9 <= MMAT( 1H 1 , 6 A6 / ) 

wkITC (o,?4> A,R,C,X2 

2 FORMAT (* A = • * F 1 2 • 1 ' , * R = *,F12.K,* C = *,F12.10, 

1* 1/(2 OMfOA SIN PHI OIST) = *,1PF12.4) 

L 

C DRAW VtLOf I TY PRCF I LES 

C 

DO 2 1L 1=1, NN 

^ 1 sox( i »=-sd( i ) 

CALL ORmW( NN,VE LO ,SDX ,1 ,7 ,L4, ITL.25. , 1( J. , 1C , 4 , 2 , 2 , 8 , l’> , 1 , L LL 1 

CALL .)PAn( NN , VELP ,SOX ,3 ,0 , L4 f ITL ,?5. , 1C 0. , 1C , 4, 2,2, tt, 1* , 1,L LL ) 

GO TC 44 
3/ CONTINUE 

S TCP 
END 
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Program 4. A program to compute the geostrophic 
surface velocity over the FNWC 63 x 63 
fields on the IBM/360 computer 
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L DIMENSION 3(6J,63),VCX{62,62)»VCY(b2,fc2) 

ill WL NS I J\ T(63,o3),VAX(62,62),VAY(62,62) ? VdX( 62 ,r^2 ) , V0Y(62,62) 
ni Vt NS ION VX (6 2,62 ) f VY( 6 2,62) , IT( 24) »D( 6 ) , I X< 62) , I Y( 62 ) 
DIMENSION V ( 6 2 , 6 2 ) 

<♦ FCRM*T (lnl, * A = ' » 1 r E 2 U • > * f 3 = ',E20.5, f C = ',620. 5/) 

5 FCkMAT (/) 
t» FORMAT ( 13, 1 3 F i.2) 

7 FCkMAT (IX, 1319/) 
k FORMAT ( 1H 1 , 2 A A *♦/ / ) 

^ K'P M AT ( 3 X , l 3 F 9 • 2 ) 

1 i FLkMAT < 3F20. 10, 16 ) 

r 

C I CUT CCNTRcL INTEGER IN COLUMN 65 
C I OUT = 0 PRINT VELGCI T Y LISTING ONLY 

C I OUT = 1 PRINT CONTOURED MAP ONLY 

C irUT = 2 PRINT 8CTH LISTING AND CONTOURED MAP 

C 

READ (5,ll)^,ti,C,I)UT 
D( 1 ) =0. 

C( 2 ) = 1 2 1 . 9 

0 ( 3 ) = 6 1 . 0 

C ( ) = 6 1 .0 

0(5) =30.6 
C ( 6 ) = 30 .5 
C ( 7 ) = 30 • 6 
C( Q ) =30 .5 

CC A 1 1= 1 , a2 

1 Y ( I ) = I 
Ri l X ( I ) = I 

C OC 10 1=1, c3 

C DC 10 J = 1 , 6 3 

0 10 S( I , J ) = 3 5 . 

00 16 K=l,d 

REAC (3) ( (T( I ,J ), 1 = 1,63) ,J = l,o3) 

READ (8) IT 
IF (K .EC. 1. ) GO TO 13 
GC TC 23 
13 DC 22 J= 1 , o2 
CC 22 1=1,62 
VX ( I , J ) =0. 

VY ( I , J ) =0 . 

V A X ( I , J ) = T ( I , J )-T( I ♦ 1, J ) 

V A Y ( I » J ) = T ( I , J ) -T ( I,J+1) 

VBX( I,J) = (T { I , J ) ♦ T ( I ♦ 1 t j ) )-. . 5 * VAX ( I , J ) 

V8Y( I,J) = (T ( I , J ) ♦ T ( I , J ♦ 1 ) )» . 5*VAY< I ,J) 

C VCX( I,J)=S( I,J)-S(IU,J ) 

C VC Y ( I , J ) =6 ( I , J ) -S ( I,J+1) 

22 CC NT I NU E 
GC ^ C 15 

2 3 DC 2 4 J=l,62 
DC 2 1 = 1, o2 

DTa = t( 1 , j )-T( j + 1, j ) 

DT Y=T( I , J )-T< I , J+l) 

B C X = ( T ( 1, J)+T< l+l,J) )^.5*CTX 
BCY= < T < I , J ) + T ( I , J +1) )- . 5-CT Y 
C DSx=S ( 1 , Jl-S ( 1 + 1, J ) 

C CS Y = 3 < I , J J-S ( 1 , J + l ) 

DAX= < VAX < I , J ) + DT X )» . 5 
OAY= ( V A Y ( 1, J ) + 0 T Y )“ .5 
CEX= { VBX( 1, J ) +BOX H . 5 
OP. Y= ( VHY ( I , J ) ♦BOY ) * .5 
C DCX= <vcx ( I , j ) ♦osx )* . 5 

C DC Y= { VC Y ( 1 . J ) + 0 S Y ) > .5 

CALI CCRfOR (I , J , CFL ) 

VX( I » J > = V X ( I,J )♦( (CAX* A ♦ DBX* B )»0( KJ-CFL ) 

VY(I,J)=VY{ I , J )♦( {DAY*t» + CBY*B ) * t)< K ) >C F|_ > 

V ( I , J ) = SORT ( V X < I , J ) *- 2 + VY( I , J 2 ) 

V A X ( I , J J =OTX 
VA Y( I , J ) =DT Y 
VB X ( I , J ) =QDX 
VE Y ( I , J ) =BUY 
C VCX( I • J ) - DS X 

C VCY(1,J)=DSY 
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2-. CCNTINUt 
NIK = 62 
NC = 62 
AW I N = - 1 00 . U 
HNU = 10. 

A / = . U0 1 
BZ = 0.0 
I JT = 0 
I C 2 N = 1 

I F ( ICUT.E0.0. OR. HUT.E0.2) GJ TO 52 
Gc TC 14 
52 CONTINUE 

W P I T t (titl A » 3 , C 

WR I T E ( o , 6 ) IT 

WR I T E ( 6 * 7 ) ( IX ( I ) ,1 = 1, 13) 

CC 31 I=l,fe2 

wR IT 2 (6,5) ( T ( I , J) , J= 1 , 13) 

WRITE (6,6) IY ( I) , ( VX< I , J ) , J= l , 13 ) 

WRITE (6,9) ( V V ( I , J ) , J = 1 , 13) 

WRITE ( 6, '-J ) ( V ( I , J ) , J= 1, l 3) 

WRITE (6,5) 

31 CONTINUE 
WRITE (6, 6) IT 

W R I T E < 6 , 7 ) ( lx< I), 1 = 14,26) 

DC 32 1=1,62 

WRITE (6,5) ( T ( I, J), 1=14, 2c) 

WRITE (o,6) I Y ( l ) , ( V X ( I,. I ) , J= 1 4 , 2 6 ) 

WRITE (c,5) (VY( I,J),J=14,26) 

WRITE ( 6 » = ) ( V < I, J) ,J=14,2t) 

WRITE (6,5) 

32 CONTINUE 

WR I T E ( 6 , 8 ) IT 

WRI TE (6, 7 ) ( IX< I ) , I =27, 3R ) 

CC 33 1=1,82 

WP ITE (6,9) (T ( I, J) , J= 27 , 3<= ) 

WRITE (6,6) I Y < I ) , ( VX< I , J ) , J = 27, 39) 

WRITE (6,4) ( V Y ( I , J ) , J = 27, 39) 

WRITE (6,4) ( V ( I,J),J=27, 3 C ) 

WR ITE (6,6) 

33 CONTINUE 
WRITE(6,c) IT 

wR ITE (6, 7) ( IX ( I ) , I =40, 52 ) 

DC 34 1=1,62 

WP ITE (6,4) < T ( I , J ) , J = 40, 52 ) 

WRITE (6,6) IY < I ) , < VX( I , J ) , J = tO,52) 

WRITE (6,9) ( V Y ( I , J ) , J = 4U, 5 2) 

WRITE (6,9) ( V ( I , .) ) » J = 40» 52 ) 

WRITE (6,5) 

34 CONTINUE 

WR I T E ( 6 , 3 ) IT 

WP ITE (6,7) ( IX ( I ) , 1 = 53, 62 ) 

CC 35 1=1,62 

WRITE (6,9) ( T ( I , J ) , J = 53 , 62 ) 

WRITE (6,6) I Y ( I ) , ( V X ( I,J) »J = 53,62) 

WRITE (6,9) <VY( I ,J ) ,J=53,62) 

WRITE (6,9) <V< I, J) ,J=53, 62) 

WR ITE (6,5) 

35 CONTINUE 

14 IE < IOUT.EO. l.CR. I0UT.E0.2) GO TO 54 
GO TO 15 

54 CONTINUE 

CALL wctmAP (V, NF,NC, IT, BNr,AZ,9Z, AMIN, IJT, ICON) 

15 CONTINUE 
STOP 

END 

SUBROUTINE CCRFQP (I,J,CFL) 

R= ( 1-32 )* ( 1-32 ) ♦< J-32 )» ( J -32) 
S=(973.71-R)/(573.71+R ) 

IF (S-. 125)4,5,5 
S= . 12 5 

5 A= 1.86603/ ( 1 . ) 

0=38 1 .0/A 
T= 1.45842E-04 
CFL= l ./ <T*S*C) 

RETURN 

ENC 
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Program 5* A program to compute the geostrophic 
surface velocity over the FNWC 63 x 63 
fields on the CDC/3600 digital computer 



o n o 
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r 

r 

r 

r 



r 



r. 

c 

r. 

c 

c 

c 



PROGRAM (,pOC'.JW ( OUTPUT , TAPE 1 * T ap£? « T APE 6 = 0UTPi IT , T a»F9 > 
DIMENSION T <63,61) ,T i , T? (3969) 

S ( 6 1 * 63 ) ,Sl ( 1969 ) *S? <3969) , Si ( 1969 ) 

A < l?) ,8 ( 1?) »n( 1?) 

VAX(63,61),VAY(63,63),\/F»X(63,63).VCY(6-),63) 
VBY (61,61) * VfX (63*63) 

VX (61 ,6 3) , VY (63 ,63) , \/ ( 63 , 6 1 ) 

W I NDx (3999) f ,.i I NOy ( l9f>9 ) ,i»flNDv(3969) 
TPaNS(39«9) * STREAM ( 1969) 

PTaHU) .DATE (3) 

T X ( 6 1 ) , T Y ( 63 ), E(5>, ET (A) . F(3), r, (S) 

( T 1 <?1 ) ,T) , ( T2 , VY ) 

(Si ,S) ♦ (S3* VX) . (9?* VCY) 



DIMENSION 
D I MENS 1 ON 
DIMENSION 
COMMON 

dimension 

DIMENSION 
DIMENSION 
DIMENSION 

integer 

EOijI VALRNCF 
EODI VALF MpF 



EDMIVALFNCE (VaX,WINDX(21)) » (VAY, WINdY (21)) 



F.OUIVALFNCE (VHX,WINDV) 

EO'J I valence ( TRANS ( 21 ) ,V) , (STREAM(^l) ,vX) 
EOUl VALFNCF (DATE* A) 

£Ol) I VAl FNr.F (P * I X ) , ( C,* I x (4 ) ) 

DATA GTAH/T70770000B/ 

DATA(RTAR(I), 1=1,*)/. 5, .45,506*0./ 



A, B AND r ARE CONSTANTS, AND ARE AVFRAGF VAl NFS TOR 
The PRESENT DEPTH of CALCULATION and that OF THF PRE- 
VIOUS INTEGRATION STEP. A AND ft aRf CHOSEN aT CON- 
STANT salinity of 34,9(0/00) 

As.OOOOo) 

B=.oooonflss 

C = - ,00074 



D IS the OEPTH INCREMENT in METERS between the LFVFL 
OF PRESFMT CALCLATiON, and the PRECEDING one 
D( 1 )=n. 

D (?.) = 30.4ft 
0(3)= 30.48 
D(4) = 30.4ft 
D (5) = 30.4ft 
D (6 ) = 60.04 
0 (7) = 60.94 
D (ft) = 121 .92 
D (9) = 34. P 
P ( 1 0 ) =200 . 

D ( 1 1 ) =20 0 . 

D (12) =200. 

ET IS AN ARRAY CONTAINING THE aODPESSES OF FNWF ETEl^S 

SURFACE 

ET(1)=I0HT SEA 
1 00ET/30 ,4flM---- 
ET (2 ) = * °HTS 1 nn 

200ET/60 ,96M 

E T ( 3 ) = 1 OhtS 200 

300F T/9] ,*4M 

ET (4) = 1 OHTS 300 

400FT/121 .RM 

ET (S) =1 OHTS 400 
60 OF T/ 1 ft? . 9M — -- 
ET (6) slOwfS 6^0 



r> r> 



C --BO0FT/243.8M---.- 

FT(7)=10 h TS R<10 

c 1200fr T /3 ft 5.R M - — 

ET (8) * 1 OhT S 12 0 0 

r 

C E IS AN ARRA Y WHICH CONTAINS THE CALLING PARAMETERS 

r FOP THE SUBROUTINE RO M STF whfn reading p NWF FIELDS 

DATE ( 1 ) *0. 

E(l) sOAtf ( 1 ) 

r Ed) is defined latep in the no loop 

E(3)=0 
E ( 4 ) =3h 
t (S) sSLTAPEl 

r 

ISVO IS A VELOCITY f i E l D printing option AS FOLLOWS 
r 0 COMPLETE VELOCITY FIELD PRINTOUT AT ALL DEPTHS 

r l SlIPEACE VELOCITY F IELO WILL «E PRINTED ONLY 

r 2 NO VELOCITY FIFLo PRINTOUT WTLL RESULT 

I S VO=2 

0 ICO IS A CONTOUR PRINT SELECTOR AS FOLLOWS 

C 0 ALL levels contoured 

1 SUQFACE COMTOUR(S) ONLY 

2 NO CONTOUR(S) printout 
ICO = 2 

IF ( ISVO.EO.2 ) GO TO 26 
r IX AND TY are arrays Which SERVE To nUmbfR The Rows 

C AND COLUMNS OF THE GRID FOR PRINTING PURPOSfS 

DO M T s i .63 
I Y ( I ) = I 
41 I X ( I ) =1 

r 

C READING in OF TFMPFRaTuRE ANo SALINIT y FIELDS FOU owe; 

26 N*ll 

DO 1 9 K s 1 , N 
IF (K.GT,R) GO TO 1?7 

60 TO ( 13S, 136, 13R, 140. 142,144, 146,1 4R) K 

135 E ( 2 ) *E T ( l ) 

1351 CALL RDmsTF(TI.E) 

IF ( E ( 3 ) .NF.01GO TO 1351 
RE AD ( 2 ) SI 
GO TO 17 
C 

136 E (2) *ET (?) 

1361 CALL RDMSTF(Tl.E) 

IF (E (3) .mf.O) GO TO 1361 
RE AD ( 2 ) SI 
RF AD ( 2 ) S3 
DO 137 I*i , 3R6R 

137 SI ( I ) a ( 1R,S2*S1 ( I ) ♦5,4fl#s3 { I ) ) /25. 

GO TO 17 

C 

13« E ( 2 ) *E T ( 3 ) 

1381 CALL RDmSTF(TI.E) 

IF (E (3) .Mf.O) GO TO 1381 
RE AD ( 2 ) S2 
00 139 1*1,3969 

139 SI < I ) * <39,04*S3 ( I > *10.96*S2 ( I ) ) /50 » 



152 



GO TO 17 



153 



140 E(2)*ET(4) 

1401 CALL RPMSTF ( T 1 * E ) 

IF (E (3) ,\JF.O)GO TO HOI 
1)0 141 Ini, 3959 

141 S 1 ( I > = fft.s6*S3 ( T > *4 1 .44 # S2 ( I ) » /5n . 

00 TO 1 7 

r 

1 42 E( 2 >»ET<S) 

1421 CALL RPMSTF(Tl.E) 

IF (£. (3) . MF , 0 ) Go TO 1*2, 

RFAD(2) 53 
DO 143 1=1.3969 

I 4 3 Si (I>*<?H.08 < *S2<T>*2l.92*S3<T> > /5(l • 

GO TO 1 7 

c 

144 E( 2 )*ET( 6 ) 

1441 call RP^STF < T l »F ) 

IF (E < 3 > .NF.°> GO TO 1**1 
RFaD ( 2 ) S? 

00 145 1=1,3969 

145 Si ( I ) = ( i7. 12* S3 ( I ) *32.«ft*S2 < T > > /5f). 

GO TO 17 

r 

1** t(2)*ET«7, 

1461 CALL RI)MSTF(Tl.F) 

IF <E (3) .MF.O) HO TO l 4ft 1 
READ (?) T? 

RFA0<2) S3 
DO 147 1=1.3959 

147 SI ( I ) * ( 156. H*S2 (I) ♦*3.84*s3 ( T ) ) /?00. 

GO TO 17 
r 

H* E (2) *ET (ft) 

1 4 ft 1 CALL RPMSTF (T1 tF) 

IF (E (3) ,NF.0)GO TO 1 48 1 
00 149 1=1,3969 

149 51 (I)=(34.24*S?(!)*liS5.7G*S3IT>>/?00. 

GO TO 1 7 
r 

127 «EAD(2) (T, (I) ,I«21,39fl9) 

RF AO (2 > Si 

r T HE FIELDS ARE WRITTFN ON DISK 

17 W0ITEH1T.S 

IF (K.EO.H) WRITE <9) t2»S8 
19 CONTINUE 
REWIND 9 
r 

c he actual calculation begins here, the first calcu- 

c LaTION INITTALI7E5 COMPONENTS 4 T THE LEVEL OF aSSUm- 

r ED HOT 7 ON 

RF AD ( 9 ) T,S 

13 00 22 jsl.52 
00 22 1=1.52 

VX < I * J> a n . 

VY(I»J)=ft. 

V ( I * J) =0 . 



.“1 o 
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VAx(I*J)aT(ItJ)-T(T*liJ) 

VAY(It J> =T < I * J> -T < T • JM > 

VRX(I*J)a(T(TtJ)*T(I*l»J)>*.5*VAX(l»J> 

VRY(I»J)a(T(I.j)*T(ItJ*l))*.5*vAY(I*j) 

VCX ( I , J) = S ( I • J) -S ( 1*1 • J) 

22 VCY(I»J>*S(I»J)-S(T»J*1> 

r 

C. THE ACCUMULATION OF COMPONENTS RFRINS HERE 
M2aN* 1 

00 15 K*?,M2 
HEAD (9) T * S 

23 uo 24 Ja 1 , 62 
00 24 I a 1 ,62 

L)TX = T < I. J) -T ( 1*1 *.i) 

0TY*T(I,j)-T(ItJ*1 ) 

HOxa ( T ( I * l)*T(I*l*J) )*.5*DTX 
BOY=<T< I* j) *T< I* J*l> >*.5 # DTY 
USX = S ( I ..)) -S( 1*1 , J) 

OSY = $ ( I* J) -S< I , J*1 ) 

DAX* (VAX (I i J) *DTX) *.5 
DAy= ( VA y ( T* J) *DT Y ) * ,5 
DR X a ( v8 X ( T . J) ♦ R D X ) »*5 
OflY* ( V0Y ( T * J) *RDy) *.5 
OCXs (VCX ( I , J) *()SX ) *.5 
t)CY= < VCY < T * J> *DSY> * .5 

C CORFOR is A SUBROUTINE FOP COMPUTING CORIOLIS FORCE 

C RASED ON GRID POSITION 

CALL CORFOR (I.J.CFL) 

VX (ItJ)xVX ( I » J) ♦ ( (DAX*A*DRX*9*0CX*C> *D ( K ) *CF|_ ) * , 4665 
VY ( I t J) aVY ( I » J) ♦ ( ( DA y # A*08Y*R*DCy«C) »D ( K) *CFL) #.4665 
1**4555 

V ( I . J) a SORT ( VX ( I * J) **?*VY ( I , J) **?> 

VAX ( I , J) xOTX 
VAy ( I , J) =OTY 
VRX < I * J) aOOX 
VRY ( I . J) =ROY 
VCX ( I « J) aOSX 
VCY ( I * J) rOSY 

24 CONTINUF 



C FOR COMP aTIBILITY ^ITH PLOTTING PROGRAM, VX, VY AND V 

C MUST BF r X TENDED TO 63*63 ARRAYS WITH DUMMY VALUES 

DO 40 1*1,53 
VX (63, I ) ryX (62,1) 

VX ( I ,63) a VX (1,62) 

VY (63, I ) rVY (62, I ) 

VY ( I , 63 ) av Y ( 1,62) 

V (63, 1 ) , v (62, 1 ) 

40 V ( I , 63 ) a v ( I # 62 ) 
o 

C CONTOUR AND PRINTING SFLEC T I ON 

10 IF ( IC0.F0.21G0 TO 16 

IF(IC0,E0.1.AND.K.NE,N)G0 TO 16 
WRITE (6,1 ) 

CALL CONTOUR (VX,63,63,3.,3lAvX) 

CALL COntOUR(VY,63,63*3.,3LAvY) 

16 IF ( ISV0.F0.2) GO TO 15 



IF ( ISVO.FO. 1 . AND.K.NF.N) GO TO 15 



155 



r 



c 



1 FORMAT(lHl) 
h FORMAT (/) 

6 FORMAT (TJ.13F9.4) 

7 FORMAT (IX, 1.319/) 

H FORMAT (3X»13F9.4) 
q format (3x.13F9,2) 



31 



WRITE <6. ] ) 

WRITE (6,7) (TX(T) *1=1*13) 

DO 31 1=1,63 

(6.9) ( T ( I » J) » J*1 • 1 3 ' 

(S ( I , J) * J*1 * 13) 

TY ( I ) , ( VX ( T * J) « J: 
(VY(I,J),«J=1,13) 
(V(I,J)»J*i»13) 
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(6,9) 
( 6 , 6 ) 
(6, A) 
( 6 , 6 ) 
( 6 , 6 ) 



WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE (6« 1 ) 
wRITE(6,7) 
DO 32 
WRITE 
WRITE 

w ritf 

WRITE 
WRITF 
WRITE 
WRITE (6, l ) 

WRITE (6, n 



1,13) 



33 



34 



(IX(T), 1=14,26) 

1=1 .63 

(6*9) ( T ( I . J) * J= l 4 . 26 > 

( S ( T * J) . j=14,26) 

I Y ( T ) , (VX(I'J) .J=]A»26) 
( VY ( I . .J) . J=14»26) 

( V ( I , J) , J* 14,26 ) 



(6,9) 
( 6 , 6 ) 
( 6 » fl > 
(6. A) 
( 6 . 6 ) 



( IX ( I ) ,1*27.39) 



1*1 ,63 



(6,9) 

(6,9) 

( 6 , 6 ) 

( 6 , 6 ) 

( 6 , 6 ) 

(6,5) 



DO 33 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE (6, 1 ) 
WRITE (6,7) 
DO 34 i *1, 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE (6*1) 
WRITE (6,7) 
DO 15 
WRITE 
WRITF 
WRITE 
WRITE 
WRITE 
WRITE 



(T ( I , J) , ja’7,39) 

(S ( I . J) ♦ J*?7,39) 

I Y ( I ) , ( V X ( T ♦ J ) * J*?7 , 39 ) 
( VY ( I , J) , J=27»39) 
<V(I,J> ,J*27,39) 



(6,9) 

(6.9> 

( 6 . 6 ) 

( 6 , 6 ) 

( 6 , 6 ) 

( 6 . 6 > 



( IX ( I ) , 1*40,52) 

63 

( T ( I , J ) ,J=40,52) 

( S ( I , J) . J=40»52> 

X Y ( I ) , ( V X ( T * J > » J = 40 ,52) 
(VY ( I , J) , J = 4 0,52) 

( V( T , J) , J=40,52) 



( IX ( T) . 1*53,63) 

1*1,63 

(6.9) (T ( I. J) ♦ J*«3»63> 

( S ( T » J) » J=63,63 ) 

I Y ( 1 ) * ( V X ( T » J ) »J = 53»63) 
( VY ( T , J) , J=53»63> 

( V ( T . J) , J*6-3 ,63 ) 



15 



(6,9) 

( 6 , 6 ) 

( 6 . 6 ) 

( 6 , 6 ) 

(6,5) 
CONTINUE 
REWIND 1 tPEW I NO 



REWIND 9 
PAUSE 1 



r> n o n n non 



Ml =0 

GO TO 4“? 
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r UNLOAD TFmPEPATURF TAPF 1 AND SALINITY TAPE 2 

r mount v»tmo tape on unit i fop Reading 

r MOUNT cDPTRANS/cURPSTRM FIELD T^PE I OUTPUT) ON UntT ? 

r FOR WRITING 

r 

r. pfad in wind components 

AH G ( 1 ) =0 

G (2) =10hu W 
G (3 ) =0 
G ( A ) c3h 
G (5) aSLTAPEl 
AA CALL RDMSTF (WlNDX.G) 

IF (G (3) ,FO. 1 ) GO TO AA 
G(?)=10HV W 
A 5 CALL RDMSTF ( WlNOY. G> 

1F(G(3).F0.1) GO TO A 5 

IT is NFCFSSaRY TO COMPUTE A WIND nRlFT transport, 
SINCE ONLY the I) AND V COMPONENTS OF THE «UNn DRIFT 

ARE read IN 
DO 47 I=?l , 3RHR 

47 WINOV ( I-?0) = (SORT ( WINDY ( I ) #*P*wlNDY ( I ) **2) ) # .4Hftb 
RFW I N o I 
Ml =M1 ♦ 1 

COMBINE GFOSTROPHIC CURRENT COMPONENTS WITH WIND FIE| D 
WITH wind FIELD COMPONENTS (USING EOUI V 4 LENT ^ToRaGD 
DO 42 i = l ,*3 
DO 42 j sl f ft3 

VX ( I , J ) s V X ( I • J ) ♦ (VAX ( I , J ) ) *, AGG5 
VY(I*J»=VY(I*J>*(\/AY(I«J> ) *.A6GS 
4? V ( I *J)=V ( I « J) ♦V/HX ( I tj) 

the transport field is now To rE writtfn on tape 

A3 E ( A ) XAO 

E (5) r5LTAPE2 
TRANS (2) =1A0 r 

TRANS (2) = ( TRANS (2) .OR. DATE ( 1 ) ) 

TRANS ( 3 ) = 1 Oh 

Trans ( h) = 10 hac ur trans 

TRANS (7) xlOHCLIMATIC 
IF (Ml ,NF . 1 > TRANS (R) =1 OwG-S TROPHIC 
IF (Ml .EG. 1 ) TRAnS (R) =1 Oh TOTAL 
TRANS ( 9 ) s 1 OH NM/DAY 
TRANS (10) = 1 Oh FIG 7. 

25 CALL WRMSTF (TRANS. F (3M 
IE (E (3) .FO. 1 ) GO TO 25 

r COMPUTE STREAM FUNCTION USING SUrroUTInE PSNF AND PSFN 

CALL PSNF (PTAH.VX. VY.GTAB.ERR) 
r 

C THE STRFAM FUNCTION FlFLO IS NOW TO HE WRITTEN On TAPE 

F ( 1 ) = 0 
F ( 2 ) = 3 « 

F (3) = 5i TAPE2 

STREAN( 2 ) -iaoh 



STREAM (2) = ( STREAM ( 2) .OR.OATE ( l ) ) 
STREAM(3>ilOH 

stream <M =iohacurrstrm 

STREAM(7)=10hClIMATIC 

if (M l . nf.i ) stream (8) *1 ohg-s trophic 

IE(Ml.EO.nSTREAM(8)«lf>HTOTAL 
STREAM(R) =10h NM/OAY 
STREAm(10)=10h fig 7. 

CALL WRMSTE (STRFAM»E) 

IF (F ( 1 ) .FO. 1 ) GO TO 28 

IF (Ml .EO. 0> GO TO 48 
END FILE ? 

REWIND ? 

STOP 

END 
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One of the important advantages of the T-S Gradient method is 
the possibility of identifying the significance of the thermal and 
haline components of the geostrophic flow. Furthermore, if the 
thermal or haline components can be neglected or correlated with 
the total geostrophic flow, then the possibility exists that only one 
or the other measurement need be taken. If only the thermal struc- 
ture is known the first order thermal contribution to the geostrophic 
current can be computed: 
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Correlating the first order thermal current with the total geostrophic 
current the correlation function can then be used to determine the 
total geostrophic current from the temperature measurements alone. 

The feasibility of accomplishing this is enhanced by the de- 
velopment of the expendable bathythermograph (XBT) which allows 
a ship to steam at full speed over a region and measure the tempera- 
ture structure to 6, 000 feet. Prior to the development of the XBT, 
temperature structure was commonly measured by mechanical 
bathythermographs which were limited to 900 feet. 

If this technique were possible, one application where it would 
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be valuable would be in the International Ice Patrol Survey of the U. S. 
Coast Guard in the Grand Banks region. The purpose of these sur- 
veys is to provide information on the currents in this region for ice- 
berg drift movement during the ice season. This information is use- 
ful for safe navigation. Unfortunately, two requirements of the sur- 
vey are conflicting. First is the requirement to survey a relatively 
large area with sufficient density to delineate the flow. Second is 
the requirement to complete the survey in a short enough time to 
give a reasonable synoptic picture. At the present time the surveys 
take about two weeks during which time approximately 100 hydro- 
graphic stations are occupied (Lenczyk, 1964). If even a 50 percent 
reduction in the number of complete hydrographic stations could be 
achieved by sampling these at full speed using XBT's the survey time 
could be reduced by several days and the currents would be far more 
synoptic than the present program allows. Unfortunately, this re- 
gion is a region of rapid changes in water structure. The charts of 
dynamic topography are considered of little value after two weeks 
and another survey must be conducted. Several surveys are needed 
in a single ice season, March through June. 

Three water masses have been identified in the Grand Banks 
region (Kollmeyer, 1966): Gulf Stream Water, Labrador Water, 
and mixed water, representing the boundary between the two other 
Waters. If the geostrophic currents are to be determined from the 
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thermal component alone some technique must be established for 
separating the water masses as a different correlation will exist be- 
tween the thermal component and the total flow in each water mass. 

For the purpose of testing the feasibility of establishing such a 
relationship 33 hydrographic station pairs in the Gulf Stream Water 
were selected from the U. S. Coast Guard Ice Patrol Reports 
(Moynihan, 1968). The stations were selected using the Coast Guard 
charts of dynamic topography for the survey containing each station 
pair. The station pairs were used to compute the thermal component 
of the geostrophic current using Equation A. Ill- 1 , and the total geo- 
strophic current. The coefficients used in Equation A. IV - 1 were: 

K = 8 3. 4 x 10' 6 
K = 8 5. 8 x 10" 7 

The resulting thermal component is plotted against the total geo- 
strophic surface velocity in Figure A. IV-1. A linear regression 
analysis fit in the least square sense, yields the following expression 

(V - V ) = 0. 39 (V - V ) + 2. 84 (cm/sec) 

I L+ 1 ^ rp 

The variance of the total geostrophic surface velocity is 15. 5 

/ 2 
(cm/ sec) . 

Therefore, if the water mass boundaries can be identified it 
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Figure AIV-1. Relationship between the total geostrophic surface current and the T- S Gradient 
thermal component in the Gulf Stream Water off the Grand Banks. 
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is possible to compute the total geostrophic surface current from the 
thermal component using the above expression. Hydrographic sur- 
veys of this region could be made more synoptic using XBT's to 
substitute for some of the standard hydrographic stations. 



164 



APPENDIX V 



Field Evaluation 



165 



FIELD EVALUATION 



Introduction 



Based on the results of this research the T-S Gradient method 
should allow the computation of surface geostrophic currents in good 
agreement with those computed using the standard geostrophic calcu- 
lations. Of course, both computations suffer from the limitations 
inherent in the geostrophic assumptions. Therefore, the purpose 
of this research to find a simplified scheme to make geostrophic 
computations has been met, yet the computation of currents without 
environmental verification is not very satisfying. However, it must 
be concluded that at this time current observations are too few and 
unreliable to be used as a quantitative measure of indirectly com- 
puted currents. This appendix reviews two attempts by the author 
to verify the representativeness of geostrophic currents to the actual 
ocean currents, and proposes some refinements of these attempts 
to continue the research. These two experiments are the only ones 
known to the author as direct efforts to verify computed surface 
currents in the open ocean in the slow California Current. 

The major oceanic gyres appear to be revealed by geostrophic 
computations (Reid, 1961, and Stommel, 1965). Yet at any point 
or time in the ocean the geostrophic current may or may not be 



representative of the instantaneous flow. As mentioned earlier, 
many forces act to drive the waters of the ocean, and generally the 
flow in the ocean will consist of both baroclinic and barotrophic 
components (Foffonoff, 1962). Therefore, the determination of the 
geostrophic current at any given time from a pair of hydrographic 
stations does not yield the existing current. However, if the hydro- 
graphic stations are sampled in a time-series fashion over a suffi- 
cient length of time to remove the non-geos trophic periodic com- 
ponents, then the resulting computed current can be expected to 
come into closer agreement with the average current over the period 
of the extended sampling. 

Defant (1950) has shown that the internal oscillations in the 
water structure can lead to significant errors in the geostrophic 
currents. Particularly important in dynamic computations is the 
change in the water structure associated with internal waves, and 
the most significant periods are those of the diurnal and semidiurnal 
tide (Defant, 1950). LaFond (1951) also identified the semidiurnal 
and diurnal periods as the most significant perterbations in the 
water structure. 

Little is known of the internal oscillations in the ocean at the 
present time. However, Haurwitz, et al. ( 1959 ) show that signifi- 
cant peaks occur in their long series of records of temperature at 
50 m and 500 m offshore from Castle Harbor, Bermuda, only at the 



semidiurnal and inertial or diurnal periods. Except for a broad 
maximum at about 0. 5 cycles per hour their spectra decreases 
monotonically with increasing frequency. There appears to be no 
statistically significant coherence between observations at 50 m and 
500 m. Cox (1962) shows in observations taken from the U. S. Navy 
Electronics Laboratory oceanographic tower with isotherm fol- 
lowers that the observations between laterally- spaced sensors are 
coherent only at low frequencies, well below the mean value of the 
Vaisala frequency. 

It would be advantageous to use more than one ship to deter- 
mine geostrophic flow between hydrographic stations, such that the 
measurements applied at the same "instant" in time. In general 
oceanographic operations are conducted by single vessels. Single 
vessel hydrographic surveys are limited to two possible sampling 
procedures, both of which may significantly influence the resulting 
computed geostrophic current. The choice lies between sampling 
between stations, first one station then another, such that enough 
casts are made at each station to remove the variability; or second 
sampling each station in a time-series fashion a sufficient number 
of times to remove the variability, moving on to the next station at 
the completion of the first. 

Rapid sampling between stations means that the samplings are 
almost synoptic, that is, taken under nearly the same conditions. 
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Furthermore, in running between stations the intervening water 
structure can be sampled. However, in the California current area 
hydrographic stations for determining geostrophic flow should be at 
least 35 kilometers apart. This means that a conventional oceano- 
graphic vessel traveling at ten knots will require two hours to 
steam between stations, and a total of four hours will be required 
to complete both stations. Changes in the hydrographic conditions 
at either station may be significant during this four hour period 
(Defant, 1950). 

On the other hand, time-series sampling at a single station 
requires a sample of sufficient length and density to remove the 
periodic variations. At least 12 hours and preferably 24 hours of 
sampling should be made at each station; with a cast every two hours, 
50 hours will be required to completely sample a station pair (24 
hours per stations with two hours steaming). One disadvantage of 
single station time-series sampling is that between any two station 
pair only one sampling of the intervening water structure is made. 

To test the two types of sampling (continuous time series 
sampling at one station versus alternating sampling between stations) 
two experiments were performed. In the first experiment the al- 
ternating sampling between stations was used. In the second, the 



24 hour sampling at each station was performed. 
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Experiment I 

The first experiment was performed from August 31, 1967 to 
September 4, 1967 off Santa Cruz, California, 

Purpose 

This experiment had a dual purpose: first, the verification of 
computed surface currents for this research; second, to examine 
the validity of Richardson's n 4/3 law 11 for particle diffusion, at large 
particle separations (in this case separations in excess of six 
nautical miles). The results of this latter work have been pre- 
viously reported (Denner, et al. , 1968). 

Sampling Procedure 

The experiment consisted of alternate sampling between two 
stations A and B at intervals of every four to six hours. The 
location of these stations were: 

A = 36* 50 1 N, 122° 36. 2 T W 
B = 36 # 20 1 N, 123° 12. 0* W 

Temperature and salinity measurements were made using 
standard hydrographic procedures. Temperatures were determined 
as the average of paired thermometers, and salinities were deter- 
mined from water samples on a Hytech model 6210 laboratory 
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salinometer. Temperature structure was measured to 1, 500 feet 
every hour on station and at three equi-distant points between stations 
with XBT casts. Surface currents were measured at Station A using 
parachute drogues (Volkmann, et al. , 1956). The drogues were 
positioned relative to a taut reference mooring, using radar. The 
initial locations of Station A and Station B, and all later locations of 
B were made using Loran C. 

Meteorological Conditions 

The meteorological conditions were favorable for the entire 
cruise. Winds were always less than Beaufort 4 and predominantly 
below Beaufort 3 from 330* true, with corresponding sea states. 

The swell were always less than five feet from about 300° true. No 
difficulty was experienced in sampling or holding station due to 
weather conditions. 

Current Measurements 

Station A was established first by setting a taut moored sur- 
face reference buoy for current measurements, diffusion studies, 
and hydrographic casts. The drogues were seeded, as nearly as 
possible, in an octagonal pattern at a range of two nautical miles 
about the reference. Also, one drogue was placed next to the refer- 
ence marker. All drogues were shallow with a 10 m line between 
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the float and parachute. The position of each drogue relative to the 
reference mooring was recorded every 30 minutes on the ships radar 
for seven hours with the exception of drogue Number 2, which was 
lost after five hours. 

The movement of all drogues was consistent over the entire 
period except for drogues 2, 4, and 8. Drogue 2 was considered 
anomalous and dropped from the analysis. The final position of 
drogue 4 and initial position of drogue 8 were anomalous and also 
were neglected. The surface current was taken as the straight line 
distance between the initial and final accepted position divided by the 
time lapse. The results of this analysis are given in Table A. V-l. 



Table A. V- 


1. Average velocity of the 


drogues in 


Experiment I. 


Drogue 


Displacement 
(n miles) 


Tim e 
(hours) 


Speed 

(knots) 


Direction 
(degrees true) 


i 


2. 0 


6. 6 


0. 28 


340 


2 


- 


- 


- 


- 


3 


1. 2 


7. 2 


0. 17 


17. 5 


4 


1 . 0 


7. 2 


0. 14 


347 


5 


1. 8 


7. 2 


0. 25 


336 


6 


1. 7 


7. 2 


0. 24 


336 


7 


2. 0 


7. 2 


0. 28 


334 


8 


1. 8 


6. 2 


0. 29 


350 


9 


1. 4 


7. 2 


0. 19 


2 
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Hydrographic Measurements 

At the conclusion of the current measurements the hydrographic 
studies were initiated with a hydrographic cast at Station A. Station 
B was established 40 nmi from station A, bearing 250® true, approxi- 
mately at right angles to the observed currents. Casts were taken 
alternately at Stations A and B. The same bottle spacings were used 
on each cast for convenience and the maximum sampling depth was 
fixed by the amount of wire out, 1, 000 m. An expendable bathy- 
thermograph was launched at each cast with the messenger, and at 
each mid-point crossing between stations. The transit time between 
stations was approximately five hours, and five casts were made at 
each station. The hydrographic studies required 60 hours. The 
hydrographic data are on file at the Naval Postgraduate School, 
Monterey, California under NPGS Cruise A-l, 1967. 

Results 



The drogue-measured current was taken as the average of the 
results summarized in Table A. V-l, 0. 25 knots toward approxi- 
mately 340° true. Geostrophic surface currents were computed 



from the hydrographic data for each successive station pair (A^B^, 
BA, A B , BA, etc. ) using both the T-S Gradient method and 

1 w w U w J 



the standard geostrophic method. A reference level of 700 m was 



Table A. V - 2. Comparison of standard and T-S gradient geostrophic surface currents, and FNWC 
currents for Experiment I. 

Standard T-S Gradient 

Station Method Method FNWC 
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selected for both calculations as this depth was achieved on each 

ca st. 

The velocity proliles indicate a steady increase in velocity 
relative to the 700 db reference level. The profiles obtained by the 
two approaches to the geostrophic computations are essentially the 
same. The relative geostrophic surface current was between 
4 cm/sec and 9. 5 cm/sec, and the values for each station pair is 
given in Table A. V - 2. These values are low compared to the meas- 
ured velocity at Station A. There are several possible explanations 
for the disagreement between the measured and computed velocity. 
First, part of the measured flow may have been in the barotrophic 
mode which is not revealed by geostrophic computations. Secondly, 
the current measurements and the hydrographic measurements were 
not concurrent, and the currents may have decreased after the 
measurements. This would have been determined by longer current 
measurements. In fact attempts were made to locate the drogues 
on each return to Station A but without success. 

Also shown in Table A. V - 2 are the surface current speed and 
direction read from the FNWC synoptic current and stream function 
charts. During this entire period the FNWC currents were to the 
south- southeast (160* to 180* true) at approximately 12. 7 cm/sec. 
This flow is in the opposite direction to the measured currents. 
However, it should be noted that the FNWC currents are the sum 
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of their mass driven and wind driven components. FNWC uses 
W it ting 1 s (1909) results (Equation 17) with a wind factor of 4. 8: 
taking the winds at 10 m/sec from 330* true, a wind drift of 15.8 
cm/sec to 330° true is computed. Therefore, it appears as if the 
FNWC current is governed mainly by the wind component which in 
this case does not lead to agreement with the measured currents. 
This is true even though the wind was persistent over the entire 
period. 

Conclusion 

This experiment shows that under these hydrographic condi- 
tions the T-S Gradient method gives essentially identical results to 
the standard geostrophic computations. These results give the geo- 
strophic current in the same direction as directly-measured cur- 
rents, but at a lower speed (4. 0 cm/sec to 9. 5 cm/sec computed 
versus 1 2. 8 cm/sec measured). Since there are many variables 
that could not be controlled, better agreement probably cannot be 
expected. 

Currents computed by FNWC for this area are in the opposite 
direction of measured flow at about 1 2. 8 cm/ sec. Since the FNWC 
current consists of both geostrophic and wind components, the wind 
component appears to be the dominant component in this case. 

Two factors might improve the confidence in these results. 
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First would be better navigation control than the Loran C provides 
m this area (+. 0, 5 nmi), and second, a longer period of current 
observations. A second experiment was designed and performed to 
correct these difficulties. 



Experiment II 

The results of the first experiment indicated that the T-S 
Gradient method of computing geostrophic component of the flow 
yielded better agreement between computed and observed currents 
than did INWC techniques. Yet the navigational control was not 
completely satisfying. On this cruise the ship’s Loran navigation 
was supplemented with Lorac navigation, a highly accurate system 
capable of providing fixes to within 0. 1 nmi. Experiment II was per- 
formed 50 nmi south of Point Arguello (Station A) using the U. S. 

Navy Pacific Missile Range Lorac system December 4, through 
December 7, 1967. 

Purpose 

The purpose was the same as for Experiment I. 

Sampling Procedure 

The sampling procedure was altered from Experiment I to 
provide two 24 hour stations (A and B) at which one cast was taken 
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approximately every two hours. Each station was marked with a 
taut moored surface marker. Drogue current measurements were 
made at each station. Station B was established perpendicular to 
the flow measured at Station A. A savonius rotor current meter 
was moored at Station B at a depth of 10 m. The locations of the 
stations wer e ; 

A = 33° 43. 8' N, 120° 45. 5 ! W 
B = 33° 28. 7 1 N, 121° 40. O' W 

Meteorological Conditions 



The winds were Beaufort 3 to Beaufort 5 over the entire cruise, 
with Beaufort 3 prevailing at Station A and Beaufort 5 at Station B. 

The winds were from 230° to 340° true over the entire cruise, pre- 
dominantly from 320 # to 340* true. The swell and sea were pre- 
dominantly from 320° true. 

Current Measurements 



Station A was occupied first and a reference mooring was 
established. Unfortunately the mooring did not hold and drifted 
slowly to the south- southeast. However, six drogues were seeded 
in an east-west line around the mooring. The movement of the 
drogues relative to the reference float was determined by radar. 
These measurements do not yield absolute currents since the 
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mooring was moving. However, the drogue positions were deter- 
mined frequently between casts by running alongside the drogues and 
noting their Lorac positions. Again the surface current was deter- 
mined by averaging the straight line movement between successive 
positions. The measurements were divided into three intervals. 
During the first interval drogues 1 and 6 did not move consistently 
with drogues 2, 3, 4, and 5 and were eliminated. During the second 
period all drogues moved consistently. During the third interval 
drogue 6 was eliminated. 

During the first period, 21:33 on December 4, 1967 to 14:44 
on December 5, 1967, the drogues moved at 9. 8 cm/sec toward 

127° true. During the second period, 14:40 to 17:40, December 5, 
1967, the drogues moved at 16. 2 cm/sec toward 129° true. During 
the last period. 17:28 on December 5, 1967 to 11:28 on December 6, 
1967, they moved 1 3. 2 cm/sec toward 163° true. 

At Station B another reference mooring was set, at right angles 
to the flow measured at Station A, and again six drogues were 
seeded. On this mooring a continuous recording savonious current 
meter was placed at 10 m below the surface on the mooring. Several 
problems were encountered at this station. Radar contact with the 
drogues was lost after only a few hours, too short a time to deter- 
mine significant movement. Attempts to position the drogues by 
running alongside and taking their Lorac position failed when the 
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Lorac signal was lost. Reliable current measurements could not be 
determined from the drogues. The current meter on the mooring 
should have provided independent measure of the currents, but the 
recorder paper failed to advance and the record was of no value. 

Hydrographic Measurements 

The hydrographic stations were established prior to the cruise. 
Hydrographic casts were to be initiated every two hours, however, 
delays of up to one hour were experienced. Only 10 casts of the 12 
planned at each station were obtained. The hydrographic station 
data is on file at the Naval Postgraduate School, Monterey, 

California, under NPGS Cruise A- 2, 1967. Again, the stations 
were approximately perpendicular to the measured drift of the 
drogues at Station A. The section established will give the geo- 
strophic current components along 160/340° true. The hydrographic 
data at each station was averaged to remove insofar as possible any 
short period variations, or measurement errors. 

In a second analysis of the geostrophic currents, the XBT 
temperatures were substituted for the reversing thermometer tem- 
peratures and the computations repeated. The purpose of this 
analysis was to see if the current structure derived from the XBT 
data would agree with the analysis from the reversing thermometer 



temperatures. 
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Results 



The velocity profiles for both sets of data are shown in Figure 
A. V-l. Since the values of the geostrophic current computed by 
both the T-S Gradient and the standard method agreed at all levels 
to within 0. 02 cm/sec only the T-S Gradient value is given. The 
surface current computed by the standard method, the T-S Gradient 
method, FNWC, and the drogue measured current at Station A are 
compared in Table A. V - 3. 



Table A. V - 3. Comparison of computed and measured geostrophic 
surface currents for Experiment II. 



Standard Method T-S Gradient FNWC Measured 

Reversing Reversing 

Thermometer XBT Thermometer XBT 



Cur rent 
Speed 
(cm/ sec) 


3. 3 


4. 5 


3. 3 


4. 5 


1 2. 7 to 
20. 5 


9. 7 to 
16. 8 


Current 

Direction 

(degrees 

true) 


160* 


160* 


160* 


160* 


150* to 
160* 


125* to 
165* 



Conclusions 



The surface geostrophic current through this section is 3. 3 
cm/sec by both the standard geostrophic method and the T-S Gradient 
method using reversing thermometer temperatures. Both methods 
yield 4. 5 cm/secwhen the XBT temperatures are substituted for the 




Figure AV- i . Comparison of geostrophic velocity profiles 
using reversing thermometer temperatures 
and XBT temperatures in Experiment II. 
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reversing thermometer temperatures. The FNWC computations 
yield a surface current between 12. 7 cm/sec and 20. 3 cm/sec with 
directions between 125° true and l65 # true. These values are in 
substantial agreement with the drogue measurements of 9. 7 cm/sec 
to 16. 8 cm/sec with directions between 150 # true and 165° true. 

If the wind drift computed according to FNWC is added to the geo- 
strophic surface currents the currents would be increased by 18 
cm/sec. This would lead to surface currents somewhat larger than 
those measured for the standard and T-S Gradient geostrophic com- 
putations. However, the drogue measured surface currents aver- 
aged over a depth interval of approximately eight meters around the 
ten meter depth of the center of the parachute. Therefore, one 
would expect the measured currents to be slightly lower than the 
current at the very surface. The influence of other factors has not 
been considered. 



Recommendations 

Introduction 

Clearly the results of these two experiments cannot be con- 
sidered conclusive. Much more experimental work is necessary to 
establish any relationship between actual currents in the ocean and 
those determined by indirect computation. Even though every effort 
was made to maintain the tightest possible experimental control 
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during these studies further control is necessary and possible. Two 
recent developments would have improved the possible control in 
these experiments. These developments are the commercial avail- 
ability of reliable continuous profiling salinity- temperatur e-pr es sur e 
systems (STP) (Brown, 1968), and satellite navigation systems 
(Warriner, 1958). In the following paragraphs another experiment is 
described, incorporating these developments and other improve- 
ments that occur to this researcher for the continuation of this work. 

The STP provides a continuous profile of temperature and 
salinity as a function of depth to at least classical Nansen bottle and 
reversing thermometer accuracy (± 0. 02 %o , +. 0. 02* C, ± 4 m) 
(Bissett Berman Corporation, 1968). The continuous profile of these 
variables also eliminates the need to interpolate discrete Nansen 
data and the associated interpolation error. Another advantage of 
the STP over the classical Nansen approach is speed, since a cast 
can be completed in less than half the Nansen cast time. Finally, a 
distinct advantage of the STP is the repeatability of the measures. 

The repeatability of the STP is two times better than the absolute 
measurement accuracy, that is, ±0. 01%o f ±0. 01*C, and ± 2 m. 
Taking advantage of this repeatability very precise difference meas- 
urements necessary for geostrophic computations are possible. 

Thus, using the STP the measurement and computational errors in 
geostrophic computations could be substantially reduced. 
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The U. S. Navy has used satellite navigation on its vessels 
since 1964 (Warriner, 1968). In 1967 the system was released for 
private users, and receivers are now available commercially which 
will yield position accuracy of better than 0. 1 nmi anywhere in the 
ocean (Magnavox Corporation, 1968). This represents at least a 
live-fold increase in accuracy over the best previous open ocean 
navigation system. 

Proposed Experiment 

The purpose of the experimental work discussed here is to 
establish the relationship between the currents at a point in the 
ocean and driving forces. Since currents exhibit a turbulent nature 
over a wide range of scales it is clear that we will never be able to 
specify the instantaneous velocity indirectly. However, for the 
purpose of contemporary forecasting this impossibility is not neces- 
sary. In the continuation of this research a modification of the 
experimental procedure used in this study is recommended. 

Stations 

Three stations should be established, initially in an area where 
the surface currents are not disturbed by topographic influences, 
such that the sections between the stations enclose a triangular re- 



gion. The primary reason for the triangular array is the fact that 
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the three stations are the minimum number forming a closed ele- 
ment. Continuity conside rations can be used for a closed element to 
help reduce ambiguities in the results. Furthermore, the geo- 
strophic currents perpendicular to the three sections formed by the 
sides of the triangular station pattern will allow vector resolution of 
the actual current. Station locations should be marked with suitable 
taut moorings. Station separation should be determined by the cri- 
terion suggested by Reed and Laird (1965) that the dynamic difference 
between the stations be five times the error in the dynamic computa- 
tions (see Chapter III). However, because of the potential repeata- 
bility of the STP and precision of the satellite navigation the station 
separation might be reduced from that used in Experiments I and II, 
keeping the transit time between stations to a minimum. 

Sampling 

Sampling of the following parameters should be performed at 
each station as outlined in the following discussion: 

1. Temperature and salinity should be measured from the 
surface to 1, 000 m. 

2. Surface currents and currents at 1, 000 m should be 
measured over the period of the study. 

3. Surface wind should be measured and recorded over the 



period of the study. 
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4. Sea conditions should be measured and recorded over the 
period of the study. 

5. Barometric pressure should be recorded over the period 
of the study. 

Sampling should be conducted such that significant periodic variation 
in the structure can be identified. 

Phase I - Establish Initial Conditions 

At each station the taut line moor should be established, 
drogues with radar transponders should be set at the surface and 
1, 000 m, and the best satellite position obtained. Once the station 
has been established, two STD lowerings should be made to 1, 000 
db, Nansen bottles with reversing thermometers placed at 0, 500, 
and 1, 000 db. During the period of the STD casts XBT's should be 
launched every ten minutes. This procedure should be followed at 
each station. The vessel should proceed between stations at maxi- 
mum cruising speed. Between stations XBT*s should be launched 
every one-half hour. Given a station separation of 40 km the first 
phase will require approximately 24 hours. Phase I establishes the 
stations to be used in Phase II, the initial hydrographic conditions, 
and the short term variability. 
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Phase II ~ Sampling 

The purpose of Phase II is to establish the median conditions 
over the period to provide the best estimate of the factors contribut- 
ing to the surface flow. Only the significant factors should be con- 
sidered. Returning to each station in initial order, 24 hour time- 
series sampling should be initiated to the 1, 000 db level. STD casts 
will be made every hour, with XBT probes launched every one-half 
hour. A parachute drogue will be established at the start of this 
phase at 1, 000 db and at the surface, the position of the drogues in 
the water at each station should be established every one-half hour. 
As the drogues are lost, they should be replaced. The ships posi- 
tion will be maintained at the reference position. Again the vessel 
should proceed between 24 hour stations at maximum cruising speed, 
launching XBT's every one-half hour. 

At the conclusion of the time series sampling the vessel should 
proceed to locate as many drogues as possible, recover the radar 
transponders and then the surface buoys. Approximately 96 hours 
would be required for Phase II exclusive of equipment recovery time. 

Analysis 

The currents through the array should be determined by both 
direct and indirect methods. The computed surface currents should 
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be compared to the measured surface currents. Wind data should be 
used to compute the wind component by various formulae to deter- 
mine the best expression to use to obtain agreement between meas- 
ured and computed surface currents. Significant periodic variations 
in the measured currents could be correlated with variations in the 
structure and local winds and tides to determine any relationships 
between these factors and the currents. 

Many other questions could be answered by the data collected 
in this experiment. These questions should not be ignored, but the 
goal of the experiment is to determine the relationships between ob- 
served and computed ocean surface currents. 



189 



DISTRIBUTION LIST 



No. of Copies 



Defense Documentation Center 20 

Cameron Station , Bldg . 5 
Alexandria, Virginia 22314 

Fleet Numerical Weather Center 20 

Monterey, California 

Dean of Research Administration 2 

Naval Postgraduate School 
Monterey, California 93940 

Library 2 

Naval Postgraduate School 
Monterey, California 93940 

Professor Warren W. Denner 25 

Department of Oceanography 
Naval Postgraduate School 
Monterey, California 93940 



190 



(Abstract, Continued) 

specific volume on temperature, a(T)g p , can be expressed by a quadratic 
equation and the dependence of specific volume on salinity, o:(S)^ p , can 
be expressed by a linear equation. The coefficients in the temperature de- 
pendence equation are a function of salinity and pressure, and the coef- 
ficient in the salinity dependence is a function of temperature and pressure. 
However, it is shown that the lack of experimental data and the presence of 
small errors in the data lead to inconsistencies in the values of the coef- 
ficients derived. This difficulty is corrected using Ekman's equation of 
state for sea water to generate P-V-T-S data. 

Using derived expressions for (da/dT)g p and (da/SS)^ p and equa- 
tion is derived from the Helland-Hansen relationship and is called the 
Temperature-Salinity Gradient scheme for computing geostrophic currents, or 
simply T-S Gradient scheme. This scheme gives equivalent results to the 
standard geostrophic computations, yet is computationally much simpler and 
faster than the standard scheme. Other advantages are (1) the variables 
observed in the ocean, temperature and salinity structure are used directly in 

the computations, (2) no other quantities such as <J and a need be 

o t 

computed, and (3) thermal and haline contributions to the geostrophic flow 
can be expressed independently. 

The T-S Gradient scheme is applied to existing hemispheric fields 
of T and S data. It is shown that in the expressions for (da/dT) and 
(da/9S)^. p the use of fixed coefficients over the entire Pacific Ocean 
relative to a reference level of 1000 db, introduces less than five percent 
error in the surface velocity. This application shows the feasibility of 
making geostrophic computations by the T-S gradient scheme. Unfortunately, 
it also shows the deficiencies of the available fields for hemispheric compu- 
tations by any scheme. These deficiencies are (1) the inaccuracy of the 
data, (2) the large grid spacing, and (3) the inconsistency of the temperature 
and salinity fields. These must be corrected before hemispheric geostrophic 
computations can be made. The construction of new fields to meet the re- 
quirements of geostrophic computations using available data is discussed. 

Two other aspects of current computation were investigated. First, 
using the capability of the T-S Gradient scheme to separate flow into thermal 
and haline components, the relationship between the thermal component and 
the total geostrophic surface current in the Gulf Stream water mass in the 
Grand Banks was studied. Results show that these quantities are satis- 
factorily related by a linear expression which allows the determination of 
the total geostrophic surface flow in this water mass using only temperature 
measurements. Such techniques may be useful in reducing survey time of 
the Coast Guard Ice Patrol in this region. 

A second aspect investigated is the relationship between the indirectly 
computed currents and measured currents in the California current. Two 
experiments were performed in which drogue-measured surface currents were 
compared to indirectly computed currents given by the Fleet Numerical Weather 
Central and those computed by the Helland-Hansen and T-S Gradient schemes 
using standard hydrographic data. This investigation points out the signifi- 
cant need for more work on this aspect. 
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